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Abstract 

A comprehensive, five-dimensional calculation of Higgs-boson production in gluon fu- 
sion is performed for both the minimal and the custodially protected Randall-Sundrum 
(RS) model, with Standard Model fields propagating in the bulk and the scalar sector 
confined on or near the IR brane. For the first time, an exact expression for the gg — > h 
amplitude in terms of the five-dimensional fermion propagator is derived, which includes 
the full dependence on the Higgs-boson mass. Various results in the literature are rec- 
onciled and shown to correspond to different incarnations of the RS model, in which the 
Higgs field is either localized on the IR brane or is described in terms of a narrow bulk 
state. The results in the two scenarios differ in a qualitative way: the h — > gg amplitude 
is suppressed in models where the scalar sector is localized on the IR brane, while it 
tends to be enhanced in bulk Higgs models. There is no smooth cross-over between 
the two scenarios, since the effective field-theory description breaks down in the transi- 
tion region, where contributions from higher-dimensional operators are unsuppressed. A 
detailed phenomenological analysis of Higgs production in various RS scenarios is pre- 
sented, including minimal models and extended models with a custodial symmetry, in 
which the scalar sector is either localized on or lives near the IR brane. For each scenario, 
the regions of parameter space already excluded by LHC data are derived. 



1 Introduction 



The discovery of a Higgs-like boson at the LHC [1, 2] marks the beginning of a new era in 
particle physics. The properties of the new particle appear to be close to those predicted for 
an elementary scalar with couplings as given by the Standard Model (SM). The hierarchy 
problem - the question about the ultra-violet (UV) sensitivity of the scalar sector and the 
stability of the Higgs potential under quantum fluctuations - is thus more pressing than ever. 
In extensions of the SM the scalar sector can be stabilized in various ways. The most popular 
solution to the hierarchy problem is low-scale supersymmetry, which protects the Higgs-boson 
mass by linking it to the masses of its fermionic partners. An interesting alternative is provided 
by models featuring a warped extra dimension [3], in which the SM is embedded in a compact 
extra dimension of anti-de Sitter space, while the scalar sector is localized on one of two branes 
bounding the fifth dimension. The fundamental UV cutoff of the model is the warped Planck 
scale, whose value near this "infra-red (IR) brane" lies in the TeV range. These models, 
introduced by Randall and Sundrum (RS), provide particularly attractive scenarios of TeV- 
scale new physics, since in addition to the hierarchy problem they also address the flavor 
puzzle and yield an attractive framework for understanding the hierarchies of fermion masses 
and mixing angles [4-6] and the smallness of flavor-changing neutral currents [7-13]. 

Precision measurements of the Higgs-boson couplings to SM particles, which are accessi- 
ble via studies of both the Higgs production cross sections and its decay rates into various 
final states, present unique opportunities to test the SM description of electroweak symmetry 
breaking and search for indirect hints of new physics. In the context of warped extra dimen- 
sions, Higgs physics has been studied by several authors [14-23]. The effect on the gg — > h 
amplitude caused by the heavy b' state, the SU(2)r partner of the top quark predicted in 
RS models with custodial symmetry, was investigated in [14]. Models in which the Higgs 
scalar is a pseudo Nambu-Goldstone boson, such as warped gauge-Higgs unification scenarios, 
were studied in [15, 21]. One finds that the result for the gg — > h amplitude only depends 
on the fundamental parameter v/f of these models, but that it is insensitive to the details 
about the spectrum of the Kaluza-Klein (KK) quarks. The authors of [16, 18] have studied 
the effect of KK resonances on the loop-induced hgg and /177 couplings by working out the 
corrections to the top- and bottom-quark Yukawa couplings induced by their mixing with KK 
states. In these papers no significant contributions from the heavy KK quark states propa- 
gating in the loop were observed, because the Yukawa interactions coupling the Higgs to two 
Z2-odd fermions (the second term in the last line of (8) below) were implicitly assumed to be 
zero. 1 The possibly large effect on the Higgs-boson couplings induced by the shift of the Higgs 
vacuum expectation value (vev) relative to its SM value, which can arise in RS models with 
custodial symmetry, was emphasized in [18]. The first complete calculation of the hgg and ^77 
couplings, in which both types of Yukawa interactions in (8) were included, was performed in 
[19]. In this paper both the production of Higgs bosons in the gluon fusion process as well 
as the main decay channels were studied in an extended RS model with custodial symmetry. 

1 The fact that there are two towers of KK quark states for every massive SM quark, which is deeply 
connected to the finiteness of the 5D loop amplitude [23], was overlooked in [17]. In order to obtain a finite 
sum for the infinite KK tower, the authors made the approximation m qn = X qn v/y2 with \ qn m 1 for the 
masses of the KK quarks, see eqs. (8) and (10) of their paper, which is incorrect. 
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It was observed that the dominant corrections to the hgg and h'yy couplings arise from the 
towers of KK quark states propagating in the loop, and that these effects are to a very good 
approximation independent of the masses of the corresponding SM quarks. The production 
rate was found to be suppressed in most regions of parameter space, while the branching 
fraction for the diphoton channel h — > 77 tends to be enhanced with respect to the SM. At 
about the same time, an independent analysis of the Higgs couplings to gluons and photons 
appeared [20], which reached the opposite conclusions. In a recent paper [23], it was shown 
that the discrepancy between the two sets of results can be traced back to a subtlety in the 
calculation of the loop-induced Higgs couplings to gluons and photons. In order to compute 
the relevant overlap integrals of fermion wave functions with the brane-localized Higgs field, it 
is necessary to regularize the Higgs profile in an intermediate step and give it an infinitesimal 
width 7] [24]. When the calculation of the gluon fusion amplitude is performed in a naive 
way, the limits of sending the regulator to zero (77 — > 0) and including an infinite number 
of KK modes (N — > 00) in the sum over virtual states do not commute. This ambiguity 
disappears once the loop calculation is performed in the presence of a consistent UV regula- 
tor, such as dimensional regularization with d < 4 space-time dimensions. For the case of a 
brane-localized Higgs sector, one then obtains the results of [19] no matter in which order the 
limits are taken. The same conclusion can be reached by using a hard UV momentum cutoff 
on the four- dimensional (4D) loop integral. The physical significance of the results found in 
[20] was not fully elucidated in [23], but the discussion in that paper suggests that they might 
refer to a certain limit of a model featuring a Higgs boson living in the bulk of the extra 
dimension. It was demonstrated that the gluon fusion amplitude receives an unsuppressed 
"resonance contribution" from high-mass KK states, which can resolve the wave function of 
the Higgs boson (see also [25]). This effect is absent for a brane-localized scalar sector. 

In the present paper, we shed new light on these issues by performing the calculation of 
the gg — >■ h amplitude in a novel way, as a five- dimensional (5D) loop calculation. In this 
way the very notion of KK states is avoided, and the only relevant limit to be considered is 
that of sending the regulator 77 of the Higgs profile to zero. In the context of dimensional 
regularization, we find that this limit can be taken either before or after performing the loop 
integration. In both cases we confirm the results obtained in [19, 23]. If the width of the 
Higgs profile is kept finite, in a way that will be specified more precisely below, we recover the 
findings of [20] . They correspond to a model with a narrow bulk-Higgs field, whose shape along 
the extra dimension can be resolved by the high-momentum modes of the RS model. We also 
extend the findings of these papers in many important ways. In Section 2, we define our setup 
and present some important remarks concerning the classification of the various RS model 
considered in our study. In Section 3, we derive an exact representation of the dimensionally- 
regularized gluon fusion amplitude in terms of an integral over the mixed-chirality components 
of the 5D quark propagator in the mixed momentum-position representation, including the 
contributions of the SM quarks and the full dependence on the Higgs-boson mass. To the best 
of our knowledge, such a result has not been presented before. Our expression holds for an 
arbitrary Higgs profile. The calculation of the 5D propagator for the case of a very narrow 
Higgs profile localized near the IR brane is performed in Section 4, with technical details 
relegated to Appendix A. In Section 5, we use these results to evaluate the gg — >■ h amplitude 
and show explicitly that taking the limit 77 — > commutes with the integration over the 4D 
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loop momentum. We prove a conjecture made in [23] for the analytic form of the contribution 
of the infinite tower of heavy KK quark states. We also present an alternative derivation of 
the same result by implementing the brane-localized Yukawa terms via appropriate boundary 
conditions in the field equations for the fermion mass eigenstates. In this approach, the notion 
of an infinitesimal regulator 77 does not appear. We also consider a generalization of the 
model in which two different Yukawa matrices enter in the 5D Yukawa interactions. We then 
discuss the changes that occur when the width of the Higgs profile is kept small but non-zero, 
corresponding to the case of a narrow bulk-Higgs field. While most of our discussion refers 
to the minimal RS model with the SM gauge group in the bulk, we generalize our results in 
Section 6 to an extended RS model with a custodial symmetry protecting electroweak precision 
observables [26-28]. Contrary to the minimal RS scenario, this model allows for masses of 
KK excitations that are in reach of the LHC [29-32]. Phenomenological implications of our 
findings in the context of recent LHC data are discussed in Section 7, where we study the 
corrections to the Higgs-boson production cross section in three different versions of both the 
minimal and the custodially protected RS model. We illustrate the magnitude of the effects 
as a function of the mass of the lightest KK gluon state and the scale of the 5D Yukawa 
matrices, and derive the regions in parameter space that are already excluded by recent LHC 
measurements. Our main results are summarized in the conclusions. 



2 Setup and classification of models 

Our focus in this work is on minimal RS models, in which the electroweak symmetry-breaking 
sector is localized on or near the IR brane. The extra dimension is taken to be an S 1 /Z2 
orbifold, labeled by a coordinate <fi G [— ir, 7r]. Two branes are localized on the orbifold fixed- 
points = (UV brane) and \<fr\ = ir (IR brane). The size r and curvature k of the extra 
dimension are assumed to be of Planck size, k ~ 1/r ~ Mpi. The RS metric reads [3] 

ds 2 = e~ 2 ° w 77^ dxV - r 2 dtf = ^ ^ dafda? - -L- dt 2 ^j , (1) 

where e~ a ^ with cr(0) = kr\<p\ is referred to as the warp factor. The quantity L = cr(n) = km 
measures the size of the extra dimension. In the second equation above we have introduced 
a new coordinate t = ee a ^\ where e = e~ a ^ determines the hierarchy between the Planck 
scale and the TeV scale, and Mkk = ke sets the mass scale for the low-lying KK excitations 
of the SM particles. 2 Our primary focus is on models where the scalar sector is localized on 
(or very near) the IR brane at t = 1, in contrast to more complicated models, in which the 
Higgs boson is a 5D field propagating in the extended bulk of the extra dimension [33-38] . 
While some of these extended models are rather appealing and deserve further investigation, 
also with regard to their Higgs phenomenology, we believe that the minimal models define 
important benchmark scenarios which should be explored first. This is not least because only 
in these cases analytic expressions for the production and decay amplitudes of the Higgs boson 
can be derived. A brief discussion of bulk-Higgs models can be found in Appendix B. 

2 The dimcnsionless variable t is related to the conformal coordinate z frequently used in the literature by 
the simple rescaling z = t/M^K = R' t. 
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Before presenting our results, we find it useful to make a few comments concerning our 
definition of a brane-localized Higgs sector, which is general enough to allow for a non-zero 
width of the Higgs profile, as long as it cannot be resolved by the modes of the theory and 
hence does not affect any observables. Recall that RS models are effective field theories with 
an inherent, position-dependent UV cutoff given by the warped Planck scale [39-43] 

Auv(t) ~ M P1 = Mpi ~ = ^ . (2) 

This accounts for the fact that they do not provide a description of quantum gravity The 
variation of the UV cutoff along the extra dimension is a crucial feature in order for RS 
models to provide a solution to the hierarchy problem. If the sector of electroweak symmetry 
breaking lives on or near the IR brane at t — 1, then the effective UV cutoff regularizing 
quantum corrections to the scalar sector is of order Ax e v ~ IOMkk- The little hierarchy 
problem is not addressed by RS models, because the theory must contain several KK modes 
(and hence the value of Axev must be in the multi-TeV range) in order to deserve the attribute 
"extra dimensional". As argued in [23], the scale Ax e v also provides the effective UV cutoff 
in loop graphs involving Higgs bosons. The condition that the modes in the effective theory 
cannot resolve the width of the Higgs boson can be stated as 

v\Y I 

7] <C - — — (brane-localized Higgs), (3) 

AxeV 

where \Y q \ sets the scale for the dimensionless, 5D Yukawa couplings of the model. Only if 
this condition is satisfied, the Higgs field can be regarded as being localized on the IR brane 
in the sense that any possible extension into the bulk does not give rise to observable effects. 
As shown in [23], another consequence of condition (3) is that the results for the loop-induced 
hgg and /177 couplings can be well approximated by performing truncated sums over a small 
number of KK modes, whose individual Yukawa couplings are evaluated in the limit 77—7-0. 
Relation (3) should be considered as a condition on the regulator 77 at fixed, physical UV cutoff 
Axev- For a brane-localized Higgs field one should take the limit 77 — > wherever possible, but 
the above condition states that keeping 77 finite but smaller than the bound on the right-hand 
side would not change the physics. 

A Higgs profile with a width 77 > ?j|Y^|/ATev must be regarded as a bulk field. On the one 
hand, the features of the Higgs profile can then be resolved by the high-momentum states in 
the effective theory. On the other hand, the equivalent relation Ax e v > v \Yq\/v shows that it is 
impossible in this case to take the limit 77—7-0. This would send the cutoff to infinity, and thus 
the model would no longer provide a solution to the hierarchy problem. In the general case, 
the gluon fusion amplitude in an RS model with a bulk-Higgs field depends in a complicated 
way on the shapes of the Higgs and fermion profiles along the extra dimension (see [20] for an 
approximate treatment; a more detailed analysis will be presented in [44]). However, we find 
that for a narrow Higgs profile, defined by the relation 

v\Y\ v\Y\ 

-7 — — <C 77 <C q (narrow bulk Higgs), (4) 

ATeV JWKK 

a model- independent expression can be derived, which generalizes the findings of [20]. 
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Model 


bulk Higgs 


narrow bulk Higgs 


transition region 


brane Higgs 


riiggs width. 


n — ^ 
V — 


v\Y q \ „ „ v\Y g \ 
AtcV " 71 ^ M KK 


^ ~ A TcV 


^ ^ A TcV 


r uwci CUIS. 


(M KK \n 
V ATeV / 


1 M KK \n 
v j?A Tc v / 

A/kK A/kk ^ A/kK ^ A/kK 

A Tc v f|V 9 | ^ »;A TcV ^ v\Y q \ 


/A/ KK \« 

U|y,| J 


(M KK \n 
v AtcV / 


Higgs profile 


resolved by 


resolved by 


partially resolved by 


not resolved 




all modes 


high-momentum modes 


high-mom. modes 




A(gg -> /i) 


enhanced 


enhanced 


not calculable 


suppressed 


Result 


model-dependent 


model- independent 


unreliable 


model-indep. 



Table 1: Comparison of the main features of various versions of RS model (see text for 
further explanation). The label "model-independent result" means that the corrections 
to the SM prediction for the Higgs production cross section can be calculated (to 
excellent approximation) without any reference to the Higgs and fermion bulk profiles. 



We will find that the results obtained under the two assumptions (3) and (4) are not 
connected to each other in a continuous way. Rather, one should consider the two scenarios 
as two different, distinguishable realizations of RS models. The situation resembles that 
encountered when one compares the original RS model, in which only gravity was allowed 
to propagate in the extra dimension while all SM fields were confined to the IR brane [3], 
with the more popular models in which all matter and gauge fields live in the bulk [6] . While 
the original model only addressed the hierarchy problem, the latter models are qualitatively 
different in that they also provide successful theories of flavor. 

The fact that there is no smooth interpolation between the two scenarios can be understood 
by investigating the structures of the corresponding effective theories in more detail. Table 1 
summarizes the main features of the various models as defined by the size of the width param- 
eter 7]. The second row in the table shows the scaling of power corrections, as represented by 
higher- dimensional operators in the effective Lagrangian of the RS model. Both in a generic 
bulk-Higgs model (with rj = 0(1)) and in models where the scalar sector is localized on the 
IR brane, effects of higher-dimensional operators in Higgs physics are suppressed by powers of 
the ratio M KK / y\ TeV , since as explained earlier the warped Planck scale A TeV is the natural UV 
cutoff of these theories. The situation changes if one considers bulk-Higgs models in which the 
width parameter t] is parametrically suppressed. Then the effective theory knows about an 
extra small parameter, and derivatives dt acting on the bulk scalar field can produce powers 
of I/77. As a result, there is a class of enhanced power corrections scaling like (M^/rjA^ey) 71 . 
In the transition region between the narrow bulk-Higgs and brane-localized Higgs scenarios, 
these enhanced power corrections become of 0(1) or larger, and hence the effective field-theory 
approach breaks down. Incidentally, as we discuss in more detail in Appendix B, a bulk-Higgs 
model with 77 1 is unnatural also from the point of view of the dynamical generation of 
such a narrow Higgs profile. 
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More specifically, consider first an idealized case, in which the electroweak scale v\Y q \ and 
the KK mass scale Mkk are of comparable magnitude (i.e., in which there is no little hierarchy 
problem). Then according to Table 1 the power corrections are parametrically suppressed in 
the bulk Higgs and brane Higgs cases, but they remain small (meaning Mkk/^Atcv <C 1) in 
the case of a narrow bulk Higgs. Even in that case, in the transition region rj ~ t>|Yg|/AT e v the 
power corrections become of 0(1), and no reliable calculations can be performed by retaining 
only a finite number of operators in the RS Lagrangian. In practice, due to the lack of KK 
modes below the TeV scale, there appears to be a little hierarchy between the electroweak 
scale and the KK scale, v \Y q \/MK_K ^ 0.3 or less. Then the power corrections in the transition 
region are even larger than 0(1), and also in the narrow bulk-Higgs case the lower bound 
on MKK/^Axev cannot be much smaller than 1. In view of this fact, one must consider the 
results derived in this paper for the narrow bulk-Higgs case with some caution. A more reliable 
calculation should stay in a regime where t] = 0(1) [44]. This has the disadvantage that the 
results will depend in a complicated way on the shapes of the Higgs and fermion profiles. If it 
turns out that this dependence is weak, however, then the results obtained here for the narrow 
bulk-Higgs scenario might serve as reasonable approximations. 



3 5D analysis of the gluon fusion amplitude 



We adopt the same definitions and notation as in the recent work [23], in which the gluon 
fusion process gg — » h was analyzed in the context of an effective 4D theory, where it is 
understood as a sum over the contributions from an infinite tower of KK quarks propagating 
in the loop. Our goal is to repeat the calculation using 5D quark propagators instead, for 
which we adopt the mixed momentum-position representation [39, 45-48] (with q = u,d) 



iS«(t, t'- p) = j d A x e* x (0\T(Q L (t,x) + Q R (t, x)) (Q L (t', 0) + Q R (t', 0))|0> 
A q LL (t, f; -p 2 ) i + A RL (t, t'; -p 2 )] P R + (L^R), 



(5) 



where P r> l = |(1 ±7s), and the symbol T denotes time ordering. We begin by considering the 
minimal RS model with the SM gauge group in the bulk. An extended model with a custodial 
symmetry will be discussed in Section 6. The minimal model contains an SU(2)l doublet 
quark field Q(t,x) and two SU{2)l singlet fields u(t,x) and d(t,x) in the 5D Lagrangian, 
each of which are three-component vectors in generation space. The 5D fermion states can be 
described by four- component Dirac spinors [4, 5] . We use a compact notation, where we collect 
the left- and right-handed components of the up- and down-type states into six-component 
vectors Ua = (Ua, ua) t and T>a = (D A , ^a) T with A = L,R, which are collectively referred 
to as Ql,r in the equation above. The Yukawa interactions of the Higgs boson with up- and 
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down-type quarks are then given by 3 



q=u,d 6 1 / n \ 

(6) 

= - E E din K^) <&\*) + h-C. , 

q=u,d m,n 

where the zeros in the diagonal blocks of the 6x6 Yukawa matrices are required by gauge- 
invariance. The function 8^(t — 1) denotes the normalized Higgs profile along the extra di- 
mension, which we take to be a regularized 5-function (see below). In the second step we have 
decomposed the 5D fermion spinors into 4D KK modes, 

Q A (t,x) = Y, Q^C*)?^^); A = L,R. (7) 

n 

The superscript n labels the different mass eigenstates in the 4D effective theory, such that 
n = 1, 2, 3 refer to the SM quarks, while n = 4, . . . , 9 label the six fermion modes of the first 
KK level, and so on. The functions Q^ R (t) denote the wave functions of the left- and right- 
handed components of the n th KK mass eigenstate along the extra dimension. The Yukawa 
couplings g^m are given in terms of the overlap integrals [19] 



Cn = ^ fdt Kit - 1) z4< m > {t) ( £ Y « ) u 



?>(t) 



Le 



(8) 



and likewise in the down-type quark sector. In the last step we have rewritten the answer 
in terms of the Z 2 -even and Z 2 -odd fermion profiles C7„ (t) and S^ft) introduced in [10], 
which are diagonal 3x3 matrices in generation space. These can be expressed in terms of 
combinations of Bessel functions, whose rank depends on the bulk mass parameters cq = 
Mq/k and c u ^ = —M u ^/k of the 5D fermion fields [4, 5]. Without loss of generality, we 
work in a basis where the c, matrices are diagonal. The SU(2)l gauge symmetry in the bulk 
implies that the 5't/(2)-doublet quark fields have common Cq parameters. The 3-component 
vectors a\f^ , on the other hand, describe the flavor mixings of the 5D interaction eigenstates 
into the 4D mass eigenstates, which are generated by the Yukawa interactions on the IR brane. 
Because of electroweak symmetry breaking, these vectors are different for A — U,D, u, d. For 
simplicity, from now on we use the generic notation Q for U, D and q for u, d. The 3x3 
matrices Y q contain the dimensionless Yukawa couplings of the 5D theory. Contrary to the 
SM, they are assumed to have an anarchical structure, meaning that they are non-hierarchical 
matrices with 0(1) complex elements. The hierarchies of the Yukawa matrices of the SM 
quarks in the effective 4D theory are explained in terms of a geometrical realization of the 
Froggatt-Nielsen mechanism in RS models [9-11, 49]. 



3 To keep the notation transparent, we do not use boldface symbols for unit and zero matrices. 
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t-2 

Figure 1: Effective hgg couplings induced by the exchange of 5D quark states. The 
positions of the vertices along the extra dimension are denoted by t± t 2 and t. 

The one- loop graph giving rise to the gluon fusion amplitude is shown in Figure 1, where at 
each vertex an integral over the fifth coordinate t = e kr ^~^ is implied, which varies between 
e = e" fcr7r ~ 10~ 15 on the UV brane and t = 1 on the IR brane. We summarize the results of 
the calculation in terms of two coefficients C\ and C5 defined by the decomposition 

A(gg ^h) = C 1 ^(U \G% CT^gg) - C 5 ^ ( \G% GT*\gg) , (9) 

where G^ a = -|e^ Q/3 G^ (with e 0123 = -1) denotes the dual field-strength tensor. Note 
that, contrary to [23], the Wilson coefficients C\ and C5 also include the contributions of the 
SM quarks. Throughout this paper, v denotes the value of the Higgs vev in the RS model, 
which differs from the SM value v^m ~ 246 GeV by a small amount [18] (see Section 7). 

In order to perform the calculation of the gluon fusion amplitude at one-loop order con- 
sistently, it is necessary to introduce two different kinds of regulators. For a brane-localized 
scalar sector, the fermion profile functions are discontinuous on the IR brane, and hence their 
overlap integrals with a ^-function type Higgs profile are ill defined. Before computing these 
integrals, it is important to regularize the Higgs profile by giving it a small but finite width 
r) <C 1 [24]. We therefore use the notation 8^(t — 1) in (6) and (8), where the regularized 
profile has unit area and support in the interval 1 — 77 < t < 1. Many of our results will be 
independent of the shape of the Higgs profile and would remain valid for the case of a general 
bulk-Higgs field, which we briefly discuss in Appendix B. Only at the end of our analysis we 
will specialize to the case of a very narrow Higgs profile, with r\ satisfying one of the conditions 
(3) or (4). Note that we use the same Yukawa matrix Y q in the two off-diagonal blocks in (6). 
For a bulk-Higgs field, the equality of the two Yukawa matrices is a consequence of 5D Lorentz 
invariance. If the Higgs field is confined to the IR brane this argument no longer applies, and 
it would in principle be possible to allow for two different Yukawa matrices Y q c and Y q ^ in 
the two terms in the last line of (8) [20, 24]. This generalization is discussed in Appendix C, 
and the corresponding results are summarized in Section 5.7. 

Secondly, as has been emphasized in [23], it is important to introduce a consistent UV 
regulator in the calculation, even though the final answer for the gluon fusion amplitude is 
UV finite. This should not come as a surprise, as it is well known that even in the 4D case 
the introduction of a UV regulator is required in order to obtain a gauge-invariant answer. To 
see this, consider the loop diagram for a single KK mode, which naively is linearly divergent. 
Using invariance under p — > —p, a superficial logarithmic divergence remains. In dimensional 
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regularization, one encounters the integral 



d d p 

(27T)° 



4- d 



P 



A 



d (p 2 -Ay (p 2 -A)\ 



e(h) ■e{k 2 ) 



(10) 



which identically vanishes for d ^ 4. Here A = m^ n — — y)m 2 h arises after combining 
denominator using Feynman parameters. Note that if the calculation was performed naively 
in four dimensions, then only the second term would be present, and it would correspond to 
a gauge-dependent operator A°\A^ ,a . In the 5D model, the UV regulator has the additional 
effect of regularizing the infinite sum over KK modes, which once again is superficially loga- 
rithmically divergent [23] . The relevant sum is of the form (recall that n = 4 labels the lightest 
KK excitation) 

3+67V „ / N 4-d 



N- 



lim 



E E 

n=u,d n=A 
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where m qn are the masses of the KK quarks and g^ n denote their effective 4D Yukawa couplings 
as defined in (8). For d = 4, one obtains different results depending on which of the two limits 
is evaluated first. However, in the presence of the dimensional regulator d < 4 the order of 
limits becomes irrelevant, and one obtains a unique answer for the sum, which in the limit 
d — > 4 (taken at the end of the calculation) coincides with the result found in [19]. 
With the regulators in place, the gluon fusion amplitude can be written in the form 



Agg ->h) = igl 



w - 1) 



x Tr 



q=u,d 
1 

V2 



Y n 







s q (t, t 2 - p - k 2 ) ${k 2 ) s q (t 2 , h- P ) t(h) s^h, t-p + h 



(12) 

where fcj denote the incoming momenta of the external gluons, a and b their color indices, and 
e(ki) their polarization vectors. We may now insert the decomposition of the 5D propagator 
given in (5) and try to simplify the result. This task is made complicated by the fact that the 
propagator functions are complicated functions of the 4-momentum p and the coordinates 
t, t' . In order to simplify the calculation, it is convenient to use in intermediate steps their 
representations as sums over KK modes. Using the KK decomposition (7), it is straightforward 
to show that 

1 



Al L (tA-P 2 ) = Yl 



P 



rat 



RL 



^— ' p l — mi 



(n) 



(t) Q^\t'), 



(13) 



and similarly for the other two propagator functions. With the dimensional regulator in place, 
the 4D loop integral as well as the infinite sums over KK modes converge, and therefore the KK 
representations provide exact representations of the 5D propagator functions. The integrals 
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over the coordinates t\ and ti of the two external gluons can then be performed using the 
orthonormality relations [10] 



J dt (t) (t) = 5 mn ; A = L,R. 



(14) 



After this is done, the 5D loop amplitude A in (12) is expressed as a single sum over KK 
modes, and we find that it can be reduced to integrals of the regularized Higgs profile with 
traces of the mixed-chirality components of the 5D propagator evaluated at t — t'. We define 



t + (pI) 



T-(pI) 



E 

q=u,d 



V 

7! 



T — 

g=u,d 



dt5 v h (t-l)Ti 



dt5 v h (t- l)Tr 



Y q 
Y} 



Y q 
Yj 



Ak(M;?l) 



A^(t,t;p|) 



2i 



(15) 



where p\ = —p 2 denotes the square of the Euclidean loop momentum after the Wick rotation. 
Matching the resulting expression for the amplitude A with the two-gluon matrix elements in 
(9), we obtain 



3 f 1 f 1 3 f 1 / m 2 \ 

Cx = - dx dy (1 - Axyy) I + {xyym 2 h ) = - J dz (1 - z) f(z) I + [z -^J 

° 5= L dx L dy ^ ^ xyv m ^ = L dz ^ J ~ ( z ^) ' 



(16) 



where is the Higgs-boson mass, x and y are Feynman parameters, and we abbreviate 
y = 1 — y and f(z) = arctanhi/1 — z. The quantities 

2 



/±K) 



e n EfL 2e 



r(2-e) 

~rri-e) 



2 2(l-e) 







g 

5p| 



T±(p 2 E 



m 



zO) 



/•oo ^ 

Jo dpEp E 2e dp^ T± (P^' m2 - 10 ) 



(17) 



are the dimensionally regularized loop-momentum integrals (after Wick rotation) over the 
functions T±(p E ) in (15), shifted by an amount m 2 . We work in the MS scheme with d = 4 — 2e 
space-time dimensions. In the last step we have integrated by parts, which is justified as long 
as the quantity pe dT±/dpE vanishes at pe = and at pe = oo. Our analysis in the following 
section confirms that these conditions are satisfied. 

In [23], we have also explored a more intuitive regularization scheme based on using a hard 
UV momentum cutoff on the loop integral. This can be readily implemented once we have 
the answers in the form given above. Setting e = and restricting the loop momentum to the 
range < pe < A, we obtain 



-m 



iO) -T ± (A 2 



m 



A 2 



d 
dA 2 



T ± (A 2 -m 2 ), 
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where A should be identified with the physical UV cutoff At c v of the RS model. 

The relations (16) are one of our main results. They provide exact expressions for the 
Wilson coefficients corresponding to the 5D loop integral. The trick of using the KK represen- 
tation in intermediate steps is legitimate and not different from similar techniques routinely 
used in 4D loop calculations. Note that in our analysis we have not taken the limit — > 0, 
which is often adopted in discussions of the gluon fusion amplitude and provides a good ap- 
proximation if the mass of the particle in the loop satisfies the inequality m 2 qn 3> m^/4. There 
would be no problem in using this approximation for the KK excitations, but for the light SM 
quarks (and to some extent even for the top quark) the Higgs mass must be kept in order to 
obtain a reliable result. The strategy adopted in [14-23] was to first evaluate the gluon fusion 
amplitude in the limit — > 0, then to subtract the contributions of the zero modes by hand, 
and finally to add back the contributions of the top and bottom quarks using the proper loop 
functions calculated with the physical value of the Higgs mass. Since in a 5D framework there 
is no distinction between zero modes and KK excitations, we are forced to keep the Higgs 
mass finite in order to include the SM contributions in the correct way. 

Our results (15) and (16) are valid for an arbitrary Higgs-boson profile along the extra 
dimension. As long as one succeeds in computing the mixed- chirality components of the 5D 
propagator in a generic bulk-Higgs model, one can use (16) to compute the corresponding 
effective ggh couplings. The limit of a brane-localized scalar sector corresponds to taking the 
limit rj — > in (15). The calculation of the function A^, L in that limit will be presented in 
the following section. It suffices to focus on one of the mixed-chirality components, since for 
space-like momenta the two components are related by A. q LR (t,t';pg) = [A. q RL (t', t; p%)]* ■ 



4 Calculation of the propagator functions A| L and 

We will now derive explicit expressions for the 5D fermion propagator in the mixed momentum- 
position representation (5). Previous studies of the warped-space 5D fermion propagator have 
been presented in [46-48]. We generalize these results by keeping for the first time the exact 
dependence on v 2 /M^ K and the full three-generation flavor structure (see also [50]), and by 
paying special attention to the effects of the regularized profile of the Higgs boson. 

The profiles Q^ R (t) form complete sets of functions on the interval t G [e, 1], subject to 
the orthonormality conditions (14). In our notation, the Dirac operator takes the form 

V=p , -M KKl5 ^-M KK M q (t) 1 (19) 



where 

M ^ = iK*Q + ^ m - l) W«) (20) 

is the generalized, hermitian mass matrix [23]. Here 5%(t — 1) denotes the normalized profile 
of the Higgs vev along the extra dimension. For a brane-localized scalar sector, we may 
without loss of generality assume that 6%(t — 1) = 8^(t — 1) are given by the same regularized 
5-function. For the general case of a bulk-Higgs field, the two profiles differ, but as described 
in Appendix B these differences vanish in the limit of vanishing rj. 
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Starting from the definition of the propagator in (5), it is straightforward to show that 

VS q (t,t' ]P ) = 8{t-t'), (21) 
where we have used the completeness relations 

Q [ A(t) Q i A ) \t') = 5(t-t'); A = L,R (22) 

n 

for the bulk profiles. For the various propagator functions, this generalized Dirac equation 
implies the coupled system of equations 

p 2 A q LL (t,t'; -p 2 ) - M KK + M q {t)j A q RL (t,t!; -p 2 ) = 6(t - t>) , 



A q RL (t,t'; -p 2 ) - M KK (-^ + M q (t)^ A q LL (t,t'; -p 2 ) = 



(23) 



and similarly for the other two functions. 4 Integrating these equations over an infinitesimal 
interval t G [f — 0, t' + 0] at fixed t' yields the jump conditions 

A RL (t' + 1 t';-p 2 )-A RL (t / -0A-p 2 ) = 



M KK (24) 
Al L (t> + 0, t; -p 2 ) - A q LL (t' - 0, t'; -p 2 ) = . 

We also need to specify appropriate boundary conditions on the UV and IR branes. In the 
presence of a regularized Higgs profile, they are 

(0 l)Al L (t u t';-p 2 ) = (l 0)A RL (t ll t';-p 2 )=0; for U = e, 1 . (25) 

This is nothing but the statement that the Z 2 -odd fermion profiles obey Dirichlet boundary 
conditions on the two branes. 

In order to solve the coupled equations (23), we first combine them to yield the second-order 
differential equation 



92 M 2 (t) dMq(t) v 2 

Mq[t) dt Pe 



Al L (t,t';-p 2 ) = -^5(t-t'), (26) 

KK 

where p\ = —p 2 /M^ K . We then solve this equation assuming that t ^ t', in which case 
the right-hand side vanishes. Next, we compute the function A q RL (t,t'] —p 2 ) from the second 
equation in (23). In the final step we determine the constants of integration by means of 
the jump conditions (24) and the boundary conditions (25). The solution of the second-order 
differential equation involves as integration "constants" functions Cj(t') with i = 1,...,8, 
which are 3x3 matrices in generation space and whose values can differ depending on whether 
t > t' or t < t' . In total, we thus have 16 functions Cf{t') and C<(t'). The jump conditions 



4 For p 2 = 0, we recover the first-order differential equations for the mixed-chirality components derived 
in [23], once we identify A 9 RL (t,t';0) = -A q RL (t,t') and A 9 LR (t,t';0) = -A q LR (t,t'). 
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impose eight relations among these functions, and the boundary conditions give four conditions 
each on the UV and IR branes. Solving these relations determines the coefficient functions 
uniquely. 

Up to this point our discussion is completely general and holds for an arbitrary bulk-Higgs 
field. Unfortunately, it is impossible to obtain a closed form of the solution for the general 
case of an arbitrary mass matrix M. q {t). Only for the special case where pe = a formal 
solution in terms of an ordered exponential can be given [23]. To proceed, we exploit the fact 
that the result of the calculation must be regularization independent in the limit r] — > 0. We 
therefore assume a particularly simple form of the regularized 5-function for the profile of the 
Higgs vev, for which we take a square box of width rj and height 1/rj: 

S ^(t-l)^-e(t-l + V ), with 77 (27) 
V M KK 

It then follows that for t < 1 — rj, where the Higgs profile vanishes, we have 

^)+«4rr , c „ (c : +1) )' 



while for t > 1 — rj we can approximate 
Ml{t) 



2 dM q (t) v 2 
q{) dt 2Ml K rf 



W \j_ hm 



Y)Y n J V v\Y c 



'KK 



'/I 



(29) 



The omitted terms are suppressed, relative to the leading one, by at least a factor 77. It will 
be useful to introduce the abbreviations 



for the positive, hermitian 3x3 matrices entering the leading term, which are given entirely in 
terms of the 5D anarchic Yukawa matrices. The general solution to (26) in the region t < 1 — t] 
is given in terms of modified Bessel functions I a (z). It can be written as 

A q LL (t,t'; -p 2 ) = V-ti 1 ^^ ° )( C ^ C2 ^ 

(31) 

^P-CQ+I^ )(C 5 (t>) C 6 (t>) 

I-c-ifot) \C 7 (t>) C s (t>) 



-1 2 



The general solution in the region t > 1 — rj can be expressed through hyperbolic trigonometric 
functions. It reads 

A , , = ( <x*h[S q P>(t - 1)] \( C 1 (f) C 2 (t>) 

LL\ 'i P ) I Q CO sh[S q 9^t-l)] \C 3 (t') C 4 (t') 

(32) 

/ sinh[5 g P>{t - 1)] W C 5 (f) C 6 (0 

^ smh[S q ^(t-l)} ) \C 7 (t>) C 8 (t>) 
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where the dependence on the coordinate t enters via the integral (for t > 1 — rj) 



1 1-t 



0"(t-l)= / dt'5^{t'-l) = , (33) 

Jt V 



and we have introduced the abbreviations 



S q = y/X* + v 2 P% , S q = y/X* + v 2 P% ■ (34) 

Once again the coefficient functions Ci(t') can take different values for t > t' and t < t' . 
Requiring that the propagator functions Az,z,(£, — p 2 ) and A^(t,t'; — p 2 ) are continuous at 
t = 1 — rj gives eight conditions, which allow us to relate the coefficients Ciit') to Ci(t'). 

In working out the solutions we neglect the infinitesimal regularization parameter 77 wher- 
ever possible, with two exceptions: First, like the profile of the Higgs vev itself, the 8 v (t — 1) 
functions vary rapidly over the interval 1 — 77 < t < 1, and hence 77 appears in (33) in an essen- 
tial way. Second, inside the quantities S q and S q the regulator appears in the product rjpE, 
and since in (17) we integrate over all values of the loop momentum there might in principle 
be contributions from very large momenta, for which rfp 2 E is comparable to the entries of X 2 
or larger. For the case of a brane-localized Higgs boson as defined in (3), such contributions 
are unphysical in view of the inherent UV cutoff of RS models, and we might therefore simply 
exclude them by hand. However, we find it more instructive to show their decoupling explicitly 
in the context of dimensional regularization. 

Further details of the solution for the coefficient functions are described in Appendix A. 
In the following section we report our final expressions for the quantities T±(p 2 E ) defined in 
(15). The dependence on the Euclidean 4-momentum enters our results via the quantities S q 
and S q introduced in (34) and via the ratio of certain linear combinations of modified Bessel 
functions, which we define as 

K-aKPe) — 7 / * \ t 7^~\ t /Mr /mi A — tJ,q. (doj 

I_ CA _i (ep E ) J CA+ i (pe) - J CA+ i {ep E ) I- CA -\ VPe) 

These quantities are diagonal matrices in generation space. A significant complication orig- 
inates from the fact that they do not commute with the matrices S q and S q , giving rise to 
non-trivial matrix products. It will be important for our discussion to exploit the asymptotic 
behavior of the ratio Ra for large and small values of Pe- Using the well-known properties of 
the modified Bessel functions I a (z), we find that for RepE S> 1 



Ra{pe) = 1 + ^ + Ca[ \X CA) + 0{pf) , (36) 
Pe ^P e 



up to exponentially small 0(e 2j3j5 ) terms. The asymptotic behavior for small values of Pe is 



F 2 (c A ) , Pe 



Ra(Pe) = i + 



p E 1 - 2c A 
where 



F 2 (c A ) + 



F\c A ) 
3 + 2c A 



Oipl) , (37) 



denotes the squared value of the profile of a chiral component of a SM fermion on the IR 
brane [4, 5]. 
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5 Analysis of the loop amplitude 



We now show how to calculate the loop integrals I±(m 2 ) in (17) for the cases of a brane- 
localized Higgs boson and a narrow bulk-Higgs field, as defined in (3) and (4). We perform 
the calculation in dimensional regularization, but we first motivate the results in the context 
of the more intuitive scheme in which a hard UV cutoff is used. We begin by collecting some 
general properties of the functions T±(p 2 E ) defined in (15), which are derived from the general 
solution to the differential equations discussed in the previous section and in Appendix A. 



5.1 Properties of the functions T±(p 



In the region of small momenta (\pe\ <C Mkk), the functions T±(p 2 E ) vary rapidly and in a 
way that is strongly dependent on the values of the bulk mass parameters q. This is expected, 
because in this momentum range their behavior is dominated by the contributions of the SM 
quarks. Remarkably, we find that at the special value pe = the results are given by the very 
simple expressions 



T+(0) = Tr l X i coth X i] = *o » T - (°) = 



(39) 



=u,d 



which only depend on the 5D Yukawa couplings, via the quantity X q defined in (30). In the 
neighborhood of this point the behavior is complicated and not described by a simple formula. 
For larger values of the Euclidean momentum, such that pe ^> Mkk, the function T + (p E ) 
converges toward a universal limiting value 



q=u,d 



X„ tanh 2X„ H 

q q 2p L 



Cq X q tanh 2X q 
cosh 2X n 



c q X q tanh 2 X q 
cosh 2X n 



Pe 



(M kk <Pb < v\Y q \/rj) 



(40) 



while T_{p E ) = 0(p E 2 ) falls off more rapidly. To derive this result, we have taken the limit 
r]p E — > and used the asymptotic expansion in (36). A dependence on the bulk mass param- 
eters enters only at subleading order. Interestingly, there exists a third region of extremely 
large Euclidean momentum, pe 3> v\Y q \/r), for which the behavior changes once again, and 
the function T + (p E ) tends to zero according to 



1 

VPe 



Y,KX 2 + 0(p E 2 ) = -^ + 



=u,d 



VPe 



(p E ^>v\Y q \/r)) 



(41) 



while still T_(p 2 E ) = 0(p E ). Note that in this region the loop momentum p E exceeds the 
value of the intrinsic UV cutoff of a consistent RS model with a brane-localized Higgs sector, 
because condition (3) implies AT e v *C v \Yq\/v- It can therefore only contribute if we consider 
a bulk-Higgs field as defined in (4). 

It follows from this discussion that the functions T±(p 2 E ) have all the properties required 
for the integration by parts in (17). The exact momentum dependence of these functions is 
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Pe/M kk 

Figure 2: Momentum dependence of the propagator function T + (p E ) for the case of 
one fermion generation and parameters corresponding to the top quark. The three 
curves refer to different values of the regulator r], as indicated. The vertical dashed 
line indicates the value of the UV cutoff of the RS model (for Ax e v = IOMrk)- 



rather complicated, and we refrain from giving explicit expressions for the general case. We 
will instead discuss the simpler case of a single fermion generation, which exhibits all the 
relevant features mentioned above. In this case we have obtained the analytic expression 



sinh 2S a 



EX 2 h(pE) S q sinh 2S q + k 2 (p E ) VPe cosh 2S q 
ud S q h(p E ) S q (cosh 25, - 1) + k 2 (p E ) V p E smh2S q + 2S q ' 1 ; 

where S q has been defined in (34), and 

ki(p E ) = 1 + R q (p E ) Rq(Pe) , k 2 (pE) = Rq{f>E) + Rq(Pe) ■ (43) 

The function Tl gcn (p E ) = vanishes trivially. It is a simple exercise to derive from (42) the 
various limiting behaviors shown in (39) - (41), simplified to the one-generation case. Figure 2 
shows the behavior of the result (42) for the parameter choices cq = —0.45, c q = 0.395, and 
\Y q \ = 2.3, which correspond to the physical mass m q = 172.6 GeV of the top quark. We 
set the KK scale to M KK = 2TeV, such that X q pa 0.2. The three curves correspond to 
different values of the regulator 77. The three regions of Euclidean momenta mentioned above 
Qoe/Mkk ~ 1, Pe/Mkk S> 1, and Pe/Mkk S> X q /r]) are clearly visible from the plot. The 
dark and light blue curves correspond to models for which Axev/^KK *C X q /r], and hence 
condition (3) defining a brane-localized Higgs field holds. The gray curve corresponds to the 
case of a narrow bulk Higgs, as defined in (4). 
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5.2 Analysis of the loop integrals I±(m 



,2^ 



Our final goal is to calculate the loop integrals I±(m 2 ) denned in (17) in the dimensional 
regularization scheme. For simplicity, however, we first consider the integral i+(0) at the 
special point m 2 = and work with a hard momentum cutoff A = Ax e v- F° r the case of a 
brane-localized Higgs sector, defined according to condition (3), we obtain from (18) 

J+(0)L „• =*,-*!- — ^ + -.., (44) 

+ ^ ^Ibrane Higgs u 1 2 A TeV 

with to and as defined in (39) and (40), respectively. The last term is a small threshold 
correction (suppressed by the UV cutoff, which we assume to be much larger than the KK 
mass scale), which is present in a hard-cutoff scheme but will not be visible in the dimensional 
regularization scheme discussed below. Such power-suppressed terms can be included via 
higher- dimensional operators in the effective Lagrangian of the RS model. Their suppression 
~ Mkk/AtgV is i n accordance with Table 1. 

The difference (to — £i) coincides with the expression for the quantity j]( CGHNP ) (summed 
over q — u,d) derived in [23] for the case of a brane-localized Higgs sector. It corresponds to 
the numerical result first derived in [19]. The same result would be obtained if one would take 
the limit 77 — >- before performing the integral over the loop momentum. For the opposite 
case of a narrow bulk-Higgs field, defined according to condition (4), the UV cutoff is such 
that the quantity T + (A 2 ) in (18) must be evaluated using (41), so that we obtain 

'+(0)1 h , k „. =t -^^ + ... (45) 

' y ' \ narrow bulk Higgs " 2 /?At y 

instead of (44). The two answers differ by an amount t\ given by the first term on the right- 
hand side in (40). The term to coincides with the expression for the quantity Eg ATZ ' ) (summed 
over q = u,d) derived in [23], which corresponds to the result first obtained in [20]. We 
emphasize that the threshold corrections are enhanced by a factor 1/rj in this case, which 
provides an example of the general behavior anticipated in Table 1 for the case of a narrow 
bulk-Higgs field. We will comment more on the structure of power corrections and the role of 
higher- dimensional operators in Section 5.3. 

It is instructive to reproduce the above results in the less intuitive, but more consistent 
(from a mathematical point of view) dimensional regularization scheme. We will argue that 
also in this case the limit of a brane-localized Higgs sector can be taken without encountering 
any ambiguities. In order to demonstrate this, we should perform the integrals over pz in 
(17) and then take the limit rj — > 0, and show that this yields the same answer as first setting 
77—7-0 and then integrating over the loop momentum. However, our explicit result in (42) and 
its generalization to three generations are so complicated that the dimensionally regularized 
integral cannot be evaluated in closed form. We will instead consider a toy model, which 
captures all important features of the exact result. To this end, we study the function 

Tr dcl (p 2 E ) = t °- tl 7 t2 + -fi= + t3 , (46) 
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which exhibits the same asymptotic behavior in the three regions as the exact result. Evalu- 
ating the integrals in (17) for this function, we obtain 

where ti/(2t 3 ) = 1 + 0(v 2 /M^ K ). While the first two contributions are associated with the 
scale Mkk, i-e. with low-lying KK modes, the third contribution is associated with the super- 
heavy scale Mkk A/) which for a brane-localized Higgs sector is larger than the physical UV 
cutoff of the RS model. Note that in the limit 77—7-0 this contribution tends to zero, leaving 
7™ odel (0) = (to — t\) as the final result for the integral after the UV regulator e has been 
removed, in accordance with (44). The same result is obtained if the limit 77 — > is taken 
in (46) before the integral is evaluated. The last term in (46) then reduces to a constant, 
which does not contribute to (17). In the dimensional regularization scheme, the case of a 
narrow bulk Higgs, for which the loop momenta can resolve the shape of the Higgs profile, is 
obtained by removing the UV regulator e at small but finite value of 77. In this case one finds 
7™ odcl (0) = t , in accordance with (45). 



5.3 Power corrections and higher-dimensional operators 

Let us add some comments concerning the size of generic power corrections, which can be 
described in terms of higher-dimensional operators added to the Lagrangian of the RS model 
(with unknown coefficients). For example, what should one expect for the magnitude of 
the leading power corrections to the Yukawa interactions (6) coupling the Higgs boson to 
bulk fermions? In general, higher- dimensional operators can be constructed by inserting one 
or more (covariant) derivatives acting on the fields. 5 These operators are suppressed by the 
fundamental, physical UV cutoff of RS models, which is of order the Planck scale. The leading 
such operators involving a fermion bilinear contain a single derivative, possibly accompanied 
by a factor sgn(0). We are thus led to study the object 

iDA"i a = -j-r- { e<J ^ $ 75 &<t> ) + terms containing gauge fields, (48) 



Mpi a M P1 V r r 

where 7° = {7^,275} are the 5D Dirac matrices and E^ denotes the vielbein [4, 5]. From 
now on we focus on the derivative terms only. Changing variables from to t, and using the 
definition of the warped UV cutoff in (2), we obtain 

Et iD Al a = (j? - 75 M KK d t ) + . . . . (49) 

Operators containing more than one derivative contain similar structures. For example, the 
5D d'Alembertian can be written as 

7T U s = ^rr ( n *-—^ d l =Tw7Y - M^ K - d t td t . (50) 



Mli Mp! \ 7-2 "*J A2j V (t) V t 



5 Note that the 5D Lagrangian does not contain any small mass parameters, which could be used to construct 
non-derivative operators of higher dimension. 
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Several comments are in order. First, we note that higher-derivative operators in the effec- 
tive Lagrangian are indeed suppressed by the position-dependent UV cutoff Auv(0> as stated 
in the Introduction. If we consider power corrections to couplings involving the Higgs boson 
(no matter whether the Higgs field is localized on or near the IR brane), the corresponding 
cutoff scale is Ax c v- The 4D derivatives contained in (49) and (50) will produce powers of 
external momenta or masses of the various fermion modes. The corresponding terms scale 
like (MKK/ATev) n - For models in which the Higgs field is a generic bulk scalar (with width 
rj ~ 1) or a brane-localized field, derivatives dt acting on the fields near t = 1 produce 0(1) 
factors, since the wave functions are naturally expressed in terms of the t variable, typically 
involving Bessel functions of argument x n t with x n = m n /M KK , or powers of t in the case of 
the SM fermions. (For a brane-localized Higgs field, these derivatives must be evaluated at 
t = 1~, i.e., by approaching the IR brane from the left.) Hence, the dt terms in the derivative 
operators shown above also give rise to (Mkr/Atcv)™ corrections. The situation changes if we 
consider a limit of a bulk-Higgs model in which the width rj of the Higgs profile becomes para- 
metrically suppressed. Then the Higgs profile itself, as well as the profiles of particles coupling 
to the Higgs field, change rapidly over a small interval of width 77 near the IR brane. In such a 
scenario, a derivative dt acting on the Higgs field or any field coupling to the Higgs boson picks 
up a factor I/77, and hence the corresponding power corrections scale like (MKK/V-^Tev) n ■ We 
thus confirm the scaling of power corrections anticipated in Table 1. 



5.4 Final expressions for the loop integrals 



The above discussion shows that in the presence of the UV regulator, and for a brane-localized 
Higgs boson, it is possible to take the limit 77 — > at the level of the functions T±(p 2 E ), before 
the loop integral is performed. Using the results from Appendix A, we extended the form (42) 
valid for one fermion generation to the general case of more than one fermion generations. For 
r] — t- 0, we find 



T + (P 2 E ) 




(51) 



where the quantity 



tanh X n 



2M 2 



Y q R q (p E )Yj 



tanh X q 

X, 



Rq(Pi 



(52) 



involves a non-trivial product of matrix-valued functions. Note that we have removed any 
reference to the matrices X q in the final expressions by using the identities Y q f(X q ) = 
f(X g ) Y q and f{X q ) YJ = Y q < f(X q ), which hold for an arbitrary function f{X q ) that has a 
non-singular expansion in powers of X 2 . 

We are now ready to derive the final expressions for the loop integrals in (17). The 
quantities T±(— m 2 — iO) computed using (51) replace the quantity t i n (44), (45), and (47), 
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while t\ has already been given in (40). Removing the UV regulator after the integral over 
the loop momentum has been performed, we obtain 

q=u,d (. 

I.(m 2 ) = E ^ Tr 

q=u,d 

where m 2 = m 2 + iO, and the function 

9( X l) Ibrane Higgs = X 1 taIlh X 1 ~ X « taIlh 2X 1 

obeys a non-singular series expansions in powers of X 2 . Note that due to the presence of 
strong-interaction phases arising from the analytic continuation from a Euclidean momentum 
p 2 E to — m 2 — iO, the functions I±(m 2 ) cannot simply be written in terms of the real and 
imaginary parts of a traces over matrices. If instead of the brane-localized Higgs boson we 
consider a narrow bulk-Higgs state, then the subtraction term t± is absent, see (44) and (45). 
The expressions in (53) remain valid also in this case, provided we use 

SWLrrow bulk Higgs = ^ q to6hX q . (55) 

The above equations are the main result of our paper. Up to some small corrections to be 
determined below, the first term on the right-hand side of the equation for I + (m 2 ) corresponds 
to the contribution of the infinite tower of KK quarks to the ggh amplitude. The remaining 
terms describe the contributions of the SM quarks. For the case of a brane-localized Higgs 
sector, the function g{X q ) coincides with an expression first obtained in [23] by means of a 
conjecture. In the present work we have derived this form. For the case of a narrow bulk- 
Higgs field, the expansion of g(X q ) to 0(X 2 ) reproduces the result derived in [20]. This 
demonstrates that the "brane-Higgs limit" considered in that paper really corresponds to the 
case of a narrow bulk scalar, as defined in (4). 
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(54) 



5.5 Alternative derivation of the result for a brane Higgs 

For the case of a brane-localized scalar sector, it has been shown in [23] that the fermion bulk 
profiles and the Yukawa couplings g^ n to the fermion mass eigenstates defined in (8) can also 
be derived in a different way, by solving the field equations for the fermion modes in the bulk 
and incorporating the effects of the Yukawa interactions by imposing appropriate boundary 
conditions on the IR brane. The Yukawa couplings are then derived by evaluating the fermion 
profiles in the limit t — > 1" (approached from the left), which defines their values on the IR 
brane by continuous extension. 

This method, which in [23] was established for individual fermion states, can also be 
applied to the infinite tower of KK modes, by imposing similar boundary conditions on the 
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5D propagator functions. Indeed, we find that with a brane-localized Higgs field the functions 
T±(p\) defined in (15) can also be computed as 



T + (P 



brane Higgs 



t.d 



V2 



Tr 



Y q 
Y} 



A^(l-,l-;p|) + A^(l-,l-;p|) 



(56) 



and similarly for T_(p 2 E ). The propagator functions A^ B are now computed by solving the 
coupled system of equations (23) without including the Yukawa term in the generalized mass 
matrix M. q {t) in (20). Instead, one modifies the boundary conditions on the IR brane, such 
that 
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1 



.V2M KK 
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KK 



1. Here 



~ _ tanhXg 

q ~ x n Yq 



(1- f;V) = (57) 



(58) 



are the modified Yukawa matrices introduced in [19]. The boundary conditions on the UV 
brane (at U = e) and the jump conditions (24) remain unchanged. It is a straightforward 
exercise to derive the propagator functions from these equations, and in particular to determine 
the mixed-chirality components at t = t' = 1~ . We have confirmed that inserting these results 
into (56) one reproduces the expressions given in (53). This method provides an independent 
derivation of the result for the brane-localized Higgs boson, in which the notion of a regulator 
r) never appears. 



5.6 Analysis of the zero-mode contributions 

We will now analyze the terms involving the matrices Z q in (53), which include the contribu- 
tions of the SM quarks, in more detail, using results derived in [10]. We first note that we can 
rewrite 

ZM) = 7^72- % R g (PE) V R q(Pe) > (59) 

KK 

with Y q as defined above. In terms of these quantities, the eigenvalue equation determining 
the KK masses reads 

det [1 + Z q {-mD] = . (60) 

The asymptotic expansion for Ra in (37) introduces the fermion profiles F(ca) next to the 
modified Yukawa matrices. We can then reexpress the answer in terms of the effective Yukawa 
matrices defined as [19] 

Y q eS = F(c Q )Y q F(c q ) = U q X q W^ where \ q = ^- m qfi (61) 

are diagonal, positive real matrices. The entries vn qi $ denote the zeroth-order values of the 
masses of the SM quarks. The unitary matrices U q and W q are defined by relation (61). 
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Including also the subleading terms in the expansion (37), we obtain 



Z q (p 2 E ) = F-\c Q )U q 



rn 



q,0 



Pe 



+ (S Q + m qfi S q m g o) + . . . 



UlF(c Q ) 



(62) 
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3 + 2c q 

F\c Q ) 
3 + 2c, 



U q X q 



(63) 



with x q = m 9i o/MKK are hermitian matrices giving rise to some small corrections of order 
v 2 /M^ K , which except for the two entries proportional to m^ 3 = m 2 carry an additional strong 
chiral suppression [10]. Introducing the abbreviation e q = Sq + m q S q m^l, and working to 
first order in v 2 /M^ K , we can rewrite the eigenvalue equation (60) in the form 



whereas 
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(64) 



(65) 



Only the diagonal elements of the matrices e q contribute when (64) and traces of (65) are 
evaluated to first order in v 2 /M^ K . It is then not difficult to show that the masses of the SM 
quarks are given by 
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q, 
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where the dots represent terms of order v 4 /M^ K and higher. Moreover, we find 
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(66) 



(67) 



(68) 



Note that while the parameters n q . are in general complex, the quantities e Qi are real. The 
sum in (67) extends over all six SM quarks. However, in practice the contributions of the light 
quarks can safely be neglected. For the third-generation quarks, we find that 



K t 



I -et 



3M£ K 



(69) 



-<V33 



up to chirally-suppressed terms, and a corresponding formula holds for This expression 
coincides with the result derived in [23]. Explicit formulae for the matrix elements (^,4)33 can 
also be found in this reference. 
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It is now a simple exercise to evaluate the Wilson coefficients C15 using (16). We obtain 

Ci = E Tr + ^1 + E Re ^) + ■ ■ ■ 



g=u,d 



3M£ K 



Re 



33 



TO 



33 



A(r t ) + A(r 6 ) + Trg{X u ) + Tr g(X d ) , 



(70) 



C 5 = ^ Im( N J B^) 
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Im 
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B(r t ) 



where = 4m|/m^ — zO, and the parameter integrals evaluate to [51, 52] 
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(71) 



In (70), 



we 



For the light SM quarks, these functions must be analytically continued to r < 1 

first present expressions that are exact up to small corrections of order v 4 /M^ K , represented 

by the dots, which are numerically insignificant. The leading effects, which involve traces over 

functions of Yukawa matrices and thus increase with the number of fermion generations, are 

exact to all orders in v 2 /M| K . The infinite sum over KK quark states contributes the trace 

term in the expression for C\. The second term contains the sum over the contributions of the 

SM quarks, whose Yukawa interactions are modified with respect to the SM by factors K qi . 

In the final, approximate expressions we have used the fact that all e qi parameters other 

than Et can be neglected to a very good approximation, and that for the term proportional 

to e t we can neglect the small deviation of the function A{r t ) ~ 1.03 from 1. Also, for the 

small 6-quark contribution, it is safe to neglect the small deviation of k& from 1. In this 

approximation, which is accurate to better than 1% for Mkk ^ 2 TeV, we observe that the 

Wilson coefficients C\ and C5 become independent of the bulk mass parameters q. They are 

entirely given in terms of the 5D Yukawa matrices of the RS model. In the SM, we have 
CfM = A ( T j + A ( T j and c sm = Q 



5.7 Brane-localized Higgs sector with different Yukawa matrices 

Before closing this section, we return to the generalization of the RS model with a brane- 
localized Higgs sector, in which one allows for different Yukawa matrices Y q c and Y q SJ< in the 
two terms in the last line of (8) [20, 24]. As discussed in Appendix C, we find that in this 
case all results discussed so far remain valid after we implement some simple substitutions. 
Instead of the matrices X q defined in (30) and Y q given in (58), we must use 

tanh X„ 



X < = 7^J Y ? Y *> Y < = ^T Y < C - (72) 



It follows that instead of (52) we now have 



KK X q 
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Also, the master formulae (53) must be generalized to read 
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J_(m 2 ) 




ReTr#(X, 



lmTrg(X q 




(74) 



Note that because X q is no longer a positive hermitian matrix, the various terms in these 
expressions can now have either sign. Finally, in the formulae for k% in (69) one must replace 
( Y ^Y U )J(Y U ) 33 by (Y U G YVY U C)J(Y U C) 33 . 

An interesting special case is that where Y q s = 0, meaning that the Yukawa couplings 
involving a product of two Z2-odd fields, given by the second term in the last line of (8), is 
put to zero. This choice was frequently adopted in the literature. It corresponds to taking the 
limit X q —7- in our results, in which case g(X q ) = 0, and the quantities K qi in (68) reduce to 
K qi — 1 — e qi . It follows that in this particular model one obtains 

C 1= Tre q + J2(l-e qi )A(r i ) + ... 

i 

A(r b ) + [1 - A(r t )] e t + e b , 



g=u,d 



(75) 



whereas C5 = 0. The first term in the first line is the result of the summation over the KK 
tower of quark states, while the second term gives the contributions of the SM quarks, whose 
Yukawa couplings are modified with respect to their values in the SM by factors (1 — e qi ). 
In practice, it suffices for all practical purposes to keep the terms shown in the second line. 
They agree with a corresponding result presented in [18]. For the top quark, the two effects 
cancel to a large extent, since 1 — A(r t ) ~ —0.03 is very small. For the bottom quark, the 
quantity e b is very small due to its chiral suppression. As a result, the overall modification of 
the Higgs-boson production rate is very small in this model. 



6 Extension to the RS model with custodial symmetry 

We will now present the generalization of the above results to an extended version of the RS 
model, in which large corrections to electroweak precision observables are avoided by means of 
an enlarged gauge symmetry in the bulk of the extra dimension. Electroweak precision tests 
are then no longer in conflict with having the masses of the lightest KK states lie in the range of 
a few TeV, in reach for direct production of these particles at the LHC [29-32] . Specifically, we 
consider an RS model based on the gauge symmetry SU (3)c x SU {2) L x SU (2) R x U (l)x X Plr- 
On the IR brane, the symmetry-breaking pattern SU(2)l x SU(2)r — > SU(2)v provides a cus- 
todial symmetry, which protects the T parameter from receiving excessively large contributions 
[26, 27]. This symmetry breaking is accomplished by means of a Higgs field transforming as 
a bi-doublet under the two 577(2) symmetries. The additional Plr symmetry, which inter- 
changes the two 577(2) groups, protects the left-handed Zbb coupling from receiving large 
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modifications [28]. On the UV brane, the symmetry breaking SU(2)r x U(1)x —> U(l)y gen- 
erates the SM gauge group. The symmetry breaking to £7(1)em is implemented by means of 
an interplay of the UV and IR boundary conditions. Thorough discussions of this model con- 
taining many technical details have been presented in [19, 53], and we will adopt the notations 
of the first paper throughout our analysis. 

The fermion representations we adopt are chosen such that they can be embedded into 
complete 50(5) multiplets used in the context of models with gauge- Higgs unification [31, 
32, 54]. As a consequence of the discrete Plr symmetry, which is instrumental in protecting 
the left-handed Zbb coupling [28] and its flavor-changing counterparts [11], the left-handed 
bottom quark has to be embedded in a SU(2) L x SU(2) R bi-doublet with isospin quantum 
numbers T\ = —T^ = —1/2. This fixes the quantum numbers of the other fields uniquely. In 
particular, the right-handed down- type quarks have to be embedded in an SU{2)r triplet in 
order to obtain an ?7(l)x-invariant Yukawa coupling. One arrives at the following multiplet 
structure for the quark fields with even Z 2 parity: 



Ql is a bi-doublet under SU(2)l x SU(2)r, while Tr transforms as (3, 1) © (1, 3). The fields 
with odd Z 2 parity have the opposite chirality. Their profiles are related to those of the Z 2 - 
even fields by the field equations. The inner and outer subscripts on the various fields denote 
their £7(1)em and U{l)x charges, respectively, which are connected through the relations 



The superscripts on the fields specify the type of boundary conditions they obey on the UV 
boundary. Fields with superscript (+) obey the usual mixed boundary conditions allowing for 
a light zero mode, meaning that we impose the Dirichlet condition (e) = on the profile 
functions of the corresponding Z 2 -odd fields. These zero modes correspond to the SM quarks. 6 
Fields with superscripts (— ) correspond to heavy, exotic fermions with no counterparts in the 
SM. For these states, the Dirichlet boundary condition is imposed on the Z2-even fields (this 
means imposing the conditions Cn^ ^(e) = on the profile functions) in order to avoid the 
presence of a zero mode. The UV boundary conditions for the profiles Sn^ ^ (£) and Cn ^ (£) 
are of mixed type and follow from the field equations. We do not explicitly show the boundary 
conditions on the IR brane, which in the presence of a regularized Higgs profile are of Dirichlet 
type for all fields, Sn (±) (l") = 0. 

Note that we have chosen the same SU(2)l x SU(2)r representations for all three quark 
generations, which is necessary if one wants to consistently incorporate quark mixing in the 

6 Note that the notation ul, cLl, u c r , Dr for these fields adopted here differs from the notation Ul, Dl, ur, 
cLr we used for the minimal RS model. 





(76) 



Y 



T R + Q x and Q = T[ + Y. 
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fully anarchic approach to flavor in warped extra dimensions. The chosen representations 
also play a crucial role in the suppression of flavor- changing, left-handed Z-boson couplings 
[11, 19]. Altogether, there are fifteen different quark states in the up sector and nine in the 
down sector. The boundary conditions give rise to three light modes in each sector, which are 
identified with the SM quarks. These are accompanied by KK towers consisting of groups of 
fifteen and nine modes of similar masses in the up and down sectors, respectively. In addition, 
there is a KK tower of exotic fermion states with electric charge 5/3, which exhibits nine 
excitations in each KK level. 

In order to simplify the notation as much as possible, it is convenient to introduce the 
vectors 



U = | U . I , u 



U' 



D = d, A = A, X=(*' 



(77) 

which collect the fields with same electric charges (2/3, —1/3, and 5/3). Upper-case (lower- 
case) symbols denote fields whose left-handed (right-handed) components are Z2 even. The 
corresponding matrices of bulk mass parameters are 

Cff = dia g( c Q> c q) > c 3 = C Q > C A = C Q > 

(78) 

Ca = diag(ev, c n , c T2 ) , c d - = diag(c T2 , c n ) , = diag(c Tl , c T2 ) , 

where each entry is a 3 x 3 diagonal matrix in generation space. Note that the fields U, D, 
and A are governed by the same bulk mass matrix cq, while u, d, and A are associated with 
three different mass matrices c u c, c T2 , and c T1 . The first two of them, c u c = c u and c T2 = c^, 
can be identified with the mass matrices appearing in the minimal RS model. The three new 
parameters contained in the matrix c Tl can be related to the other ones by extending the Plr 
symmetry to the part of the quark sector that mixes with the left-handed down-type zero 
modes, by requiring that the action be invariant under the exchange of the fields D' and D 
[19]. This extended version of the Plr symmetry implies 

c ri c r2 , (^9) 

and hence the number of independent bulk mass parameters is reduced to the same number 
as in the minimal RS model. Whether or not this equation holds will turn out to be largely 
irrelevant to our discussion. 

In generalization of (6), we now collect all left- and right-handed fields in the up, down, 
and exotic sectors into the 15-component vectors (Ua,ua) t and the 9-component vectors 
(Da, dA) T and (Aa, Xa) t (with A = L,R), to which we will collectively refer as Ql,r- The 
Yukawa couplings of the Higgs boson to these fields can then be written in the form 

C hm (x) = - J dt8l(t-l)h(x)Q L (t,x)^(^ Y ^Q R (t,x)+h.c., (80) 

q=u,d,\ 6 <? ' 
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where 

¥lV Y* = Y x =(Y d Y d ) (81) 

denote the corresponding 6x9 and 3x6 Yukawa matrices. Note that the 3x3 block matrices 
Y q appearing in these expressions are the same as in the minimal RS model. Even though the 
extended RS model with custodial symmetry has a much richer structure than the minimal 
model, it thus features the same number of parameters in the fermion sector, once relation 
(79) is imposed. 

With all the notation in place, we are now ready to generalize the analysis presented in 
the previous sections to the extended RS model with custodial symmetry. Since the Yukawa 
interactions (80) have the same structure as in (6), and since the boundary conditions on the 
IR brane are the same as in the minimal model, the only difference in the solution of the 
differential equations (23) concerns the UV boundary conditions imposed on the propagator 
functions. While the boundary conditions for fields with superscript (+) give rise to the par- 
ticular combination of Bessel functions defined in (35), R^\p E ) = Ra(Pe), the corresponding 
linear combination for fields with superscript (— ) is given by 



D (-),. v h A ~\(tpE) I- ca+ i(Pe) - I. CA+ i 2 (ep E ) I ca -i(Pe) 



■ (82) 



I CA -i{ep E ) I- ca .i 2 {Pe) ~ I- CA+ i(ep E ) I ca+ i($e) R { ^\p E 

Apart from this effect, we find that the central results (53) remain valid if we extend the sum 
over flavors appropriately, i.e. 




!X f ( Z q i-m 2 ) Zl(-m 2 ) 



83) 



sinh 2X 9 V 1 + Z^-m 2 ) 1 + z\{-m 2 

and similarly for J_(m 2 ). In analogy to (52), the matrices Z^(p 2 E ) are given by 

„ / 9 n v 2 tanh Xff „ , „ , „ , , ,+ tanh Xg „ . „ . .„ A , 

Zq<PE) = 2M^ ~^C~ Y * R ^ Y ? — X~ R Q {Pe) ' (84) 



where 

Rfj = diag( J R^ ) , r!q ) , R 3 = R ( q\ R^ = R Q , 

Ra = diag(J?it } , R£\ R^) , R d = diag(i4+\ R^) , R x = dmg(R^\ R^) , 

(85) 

which resembles the structure of the bulk mass matrices in (78). For simplicity of notation, 
we have omitted the argument p E of the various R^ matrices. 

The relevant squared Yukawa matrices entering the quantities Xg in (83) and (84), which 
are defined in analogy with (30), are given by the 6x6 matrix 

yytJ YJti+YiY* Y u Y^-Y d Y\ \ ( 2Y d Yl \ f 
aa {Y u Yj-Y d Yj Y u Yj+Y d Yn I 2Y u Yj 



with 




(86) 
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and the 3 x 3 matrices Y-Y1 = Y^Yl = 2Y d Y}. It follows that 



d d 



A A 



Trg(X<r) = Trg(V2X u ) +3Trg(V2X d ) 

q=u,d,\ 



(87) 



where the final answer is now expressed in terms of traces over the same 3x3 matrices X q as 
in the minimal RS model. 

Our next task is to reduce also the second term in (83) to traces over 3x3 matrices. From 
the definition (84), it is straightforward to derive that 
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where again we have omitted the argument p E of the matrices on the right-hand side 
of the equations. In the custodial model, the modified Yukawa matrices are defined a,s Y q = 
[tanh(y/2X q ) / (\^2X q )^ Y q [19] , with an extra factor of y/2 inserted compared with the minimal 
model. In (16), we need to evaluate the result (83) for values \p E \ -C M| K . Using the expansion 
in (37), we obtain after a straightforward calculation (again with x q = tji^o/Mkk) 
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(89) 

where V C km = Ul U d is the CKM mixing matrix. The terms shown explicitly above are of 
leading and subleading order in v 2 /M^ K . To this order, the quantity Z^(p 2 E ) vanishes. The 
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quantities are generalizations of the matrices 8a given in (63). They are defined as [19] 



u 



x u W, 



1 



1 2 c u V ( c u 



i + 



3 + 2c u 



+ VcKM W" 



d 2F 2 (c T2 ) [F 2 (-c 



1 



1 - 2c Q \ 
1 




1 - 2c T2 
1 / 


' 1 



F 2 (~c T , 
2c Q 



T 2 J J 



F 2 (~c Q )\ 



W d x d V^ KM , 
F 2 (cq)_\ 



1 - 2c 



T"2 



F 2 (- 



1 + 



5 + 2c Q J 



2c r2 



Uu x u i 
W d x d . 



(90) 



l-2c Q V^ 2 (cq) 
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U d x d 



After a lengthy calculation, we find that in analogy with (67) 
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(91) 



where 

«« = 1 - e« - | K ^(cq) 2XJ F-^cq) C/J (92) 

O L Jb 

now contains an extra factor of 2 in the last term compared with the result (68) for the minimal 
model, while 

% = + (*,)« • (93) 

We are now ready to present our final expressions for the Wilson coefficients C\ and C5 
in the RS model with custodial symmetry. To an excellent approximation, we obtain instead 
of (70) 



C, 



2v 2 Re ( Y uYjY u ) 



3M 2 K 



2w s 
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Im 



W33 
(YY^Y, ) 



A(r t ) + A(r 6 ) + Tr^(V2X u ) + 3Tr^( v / 2X d ) , 



w733 



(94) 



which once again is independent of the bulk mass parameters q. We find that this approxima- 
tion is accurate to better than 2% for Mkk <^ 2 TeV. Whereas the small corrections parame- 
terized by K qi and e qi have only a minor impact, the main difference between the minimal and 
the custodial RS models consists in the different multiplicity factors in the trace terms in (70) 
and (94). Since the functions g(X q ) start with a quadratic term, we must compare X 2 + Xj in 
the minimal model with the combination 2X 2 + 6Xj in the custodial model. Since we assume 
that the 5D Yukawa matrices in the up- and down-type quark sectors are random matrices 
of similar magnitude, it follows that the effect of the KK modes in the custodial model is 
approximately four times as large as in the minimal model. 7 



7 Based on a naive counting of degrees of freedom, this factor was estimated as 11/4 (instead of 4) in [22] 
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7 Phenomenological implications 



We now present a numerical study of our results for both the minimal RS model and its 
extension with custodial symmetry. In each case, we distinguish the two cases of a brane- 
localized scalar sector and a narrow bulk-Higgs scenario. At the end of this section, we also 
discuss the generalization of the brane-Higgs scenario with two different Yukawa matrices, 
which was discussed in Section 5.7. 

Based on the expressions obtained in Sections 5 and 6, we evaluate the Higgs-boson pro- 
duction cross section via gluon fusion relative to the SM cross section [23], 



„ °(gg ^)rs \K g \ 2 + \K g5 \ 2 

tih = —, r~\ — = 5 ' l ytl J 



where K g and K g $ parametrize the values of the Wilson coefficients, normalized to the SM value 
CfM = A ( T j + A ^ Th) ^ guch that Kg = Cl /Cf M and K g5 = \C h /Cf M . The quantity k v in (95) 
denotes the shift of the Higgs vev v in the RS model relative to the value v^m of the SM [18]. 
We determine k v from the shift to the Fermi constant derived in the RS model by considering 
(at tree level) the effect of the exchange of the infinite tower of KK gauge bosons on the rate 
for muon decay. 8 Using the definition v^m = ("\/2GV)~ 1//2 along with results derived in [10], 
we then obtain to first order in v 2 /M^ K 

I 1 , Lm w i -, , Lm w (Q(] \ 

^IminimalRS- ~ i + 4M 2 K ' 1 custodial RS ~ ^ ~ + 2 M^ K ' 1 J 

where mw — g 2 v 2 /4 is the lowest-order expression for the mass of the W boson, and L = 
-lne = ln(M P1 /A ToV ) ~ 33-34 for A TeV ~ 20-50TeV. In the custodial RS model with P LR 
symmetry, the correction is twice as large as in the minimal model. 

Concerning the contributions from the infinite towers of KK quarks to C\ and C5, we need 
to evaluate the traces of the functions g(X q ) defined in (54) and (55), which can be expanded 
in a power series in the positive matrix X 2 = 0(v 2 /M^ K ). Keeping only the first term in this 
expansion, one encounters the quantity 

7) 2 NS 7) 2 TV 2 ?/ 2 

TrX 2 = — VlfYV I 2 ^— 9 max (97} 

irA « 2M 2 K ^ i l ™ 2M 2 K 2 ' {97} 

where N g = 3 is the number of quark generations. Working with anarchic 5D Yukawa ma- 
trices, we assume that the entries (Y q )ij are random complex numbers, which with equal 
probability can take any value in the complex plane inside a circle of radius y max . The upper 
bound is imposed to satisfy the requirement of perturbativity in the Yukawa sector, and it is 
conventional to choose a value y maK ~ 3 [9]. In the last step we have used that on average 
(\{Yq)ij\ 2 ) = y^jl? It follows that, to good approximation, the effect of the KK tower 



8 If one uses instead the shift on the value of the VF-boson mass, one finds some additional contributions 
not enhanced by a factor of L, which are numerically insignificant [10]. 

9 In [23] the modulus and phase of the elements of the Yukawa matrices were chosen as random variables, 
in which case \{Y q ) t] \ 2 = J/^ ax /3. 
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of quark states scales proportional to the number of quark generations squared. While each 
entry of the Yukawa matrices Y q is a randomly distributed complex number, the central limit 
theorem implies that the sum over the Ng = 9 positive numbers in (97) is (approximately) 
normally distributed about the average value shown in the equation. It is this fact which 
allows us to predict the Higgs-boson production rate to good accuracy in terms of only the 
two parameters M KK and y max (see also [22, 23]). 

For the numerical analysis, we generate three sets of 5000 random and anarchic 5D Yukawa 
matrices, whose entries satisfy |(l^)jj| < y m ax with ?/ max = 0.5, 1.5, and 3. As a further 
constraint, we impose that these matrices correctly reproduce the Wolfenstein parameters p 
and r\ of the unitarity triangle (see [10] for explicit formulae). This requirement helps to 
eliminate some outliers in the plots presented below. We also require that, with appropriately 
chosen bulk mass parameters q, one can reproduce the correct values for the masses of the 
SM quarks; however, imposing this condition only has a minor impact on our results. 

Figure 3 shows the results for the ratio Rh defined in (95) in the minimal RS model for 
the scenarios with a brane- localized Higgs boson (left) and a narrow bulk-Higgs field (right), 
in dependence of the mass M (i) of the lightest KK gluon state. We use the mass of the 
first excited gluon state as a reference, because it is more physical than the KK scale Mkk, 
and because its value M g {i) « 2.45 Mkk is a model-independent prediction of the RS models 
considered in this work. The green, red, and blue scatter points refer to the three different 
values of y max - They have been obtained using the approximate expressions for the Wilson 
coefficients given in (70), but at the scale of the plots they are indistinguishable from the 
results one would obtain using the exact expressions in (16) and (53). We use = 126 GeV 
for the mass of the Higgs boson, and m t = 172.6 GeV and mj = m,b(mh) = 2.9 GeV for the 
masses of the third-generation quarks. While for the heavy top-quark it is appropriate to use 
the pole mass, a running mass should be used for the 6-quark. We observe that the ratio Rh 
is strictly below 1 for the case of a brane-localized Higgs sector, while it is larger than 1 for 
the case of a narrow bulk-Higgs state. This observation allows for a clear distinction between 
the two scenarios. Only for very small y m£LX , a few points exist for which Rh lies slightly below 
1. This effect is due to the modification of the Higgs vev in the RS model, which always gives 
rise to a negative contribution. 

In order to compare our predictions with experiment, we consider the cross section for the 
process pp — > h — > ZZ^*> — > Ai measured at the LHC. Since gg — > h is the dominant production 
channel, accounting for about 90% of the events in the SM, and because corrections to the 
hZZ coupling in RS models are in general very small [19, 22], we assume that any deviation 
of the rate for this process from its SM value can be traced back to new-physics contributions 
to the gluon fusion amplitude. The ATLAS and CMS collaborations have recently reported 
updated results for the ratio fizz = cr(pp — > h — > ZZ^)/a SM (pp — >■ h — > ZZ^*'), which 
were obtained using the full data set collected up to the end of 2012 (approximately 25 fb _1 ). 
The observed values are /4| LAS = 1.7 (at m h = 124.3 GeV) and /i§f s = 0.91 toll (at 
mh = 125.8 GeV) [55], 10 which we naively average to obtain fizz = 1-12 ifj^r The la range 
corresponding to this result is shown by the blue band in the two plots. In our analysis we 

10 At m h = 125.5 GeV, the ATLAS result is shifted to Mzp AS = L5 ± °- 4 > which is closcr to the CMS value 
and gives rise to the average result fizz = 1-09 j^ai- 
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M gW [TeV] M gW [TeV] 

Figure 3: Predictions for the ratio Rh in the minimal RS model, for the cases of a 
brane- localized Higgs boson (left) and a narrow bulk-Higgs field (right). The green, 
red, and blue scatter points correspond to model points obtained using y max = 0.5, 1.5, 
and 3, respectively. The overlaid lines show fits the various distributions. The area 
colored in blue represents the experimental la band (see text for further explanation). 



will assume that fizz ~ Rhi i-e., that any possible deviation from 1 is due to a modification of 
the production cross section of the Higgs boson in gluon fusion. Model points falling outside 
these bands are excluded at the 68% confidence level (CL). While for small y max = 0.5 most 
model points are in agreement with the data, it is interesting to observe that for larger y max 
the data already disfavor KK gluon masses in the low TeV range. The discrepancies between 
theory and experiment are stronger for the brane-Higgs model, because the mild tendency of 
an enhanced production rate seen in the data is in conflict with the suppression of the cross 
section predicted in this case. 

The overlaid, solid lines in Figure 3 show fits to the various distributions of model points. 
In regions of parameter space where the deviations of Rh from 1 are modest enough in order 
to be compatible with the data, a good approximation to these curves can be obtained by ap- 
proximating the functions g(X q ) in (54) and (55) by the first terms in their Taylor expansions 
and exploiting the anarchy of the 5D Yukawa matrices. In this way we find 



Rh 



v 2 



Z1V1 KK 



(±4iv| + ^ 9 -0(ITO.,| 2 ) + Lm - 



V 2 



(9* 



where the upper sign corresponds to the brane-localized Higgs sector and the lower sign to the 
narrow bulk-Higgs scenario. For randomly chosen complex elements of the Yukawa matrices, 
it follows that (|(^)ij| 2 ) = y 2 l3uX /2. The terms in brackets then evaluate to approximately 
[21.3 y^ ax + 3.6] for the RS model with a brane-localized Higgs, and [— 14.7y^ ax + 3.6] for 
the model with a narrow bulk Higgs (with L = 33.5). Relation (98) exhibits the quadratic 
dependency on the number of quark generations N g and on the maximum absolute value y max 
imposed on the entries of the random Yukawa matrices. 
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Figure 4: Excluded regions of parameter space in the minimal RS model, for the cases 
of a brane- localized Higgs boson (left) and a narrow bulk-Higgs field (right). The 
vertical dashed line shows the lower bound on M g {r, obtained from a tree-level analysis 
of electroweak precision observables. 



Even at the present level of precision, the existing measurements of the Higgs-boson pro- 
duction cross section already provide highly non-trivial constraints on the parameter space 
of RS models. In Figure 4, we show the regions in the M ff (i)-y max parameter space which 
are already excluded by the current experimental data at various confidence levels. To obtain 
these regions, we first fit an approximately gaussian distribution to the model points shown 
in Figure 3 for each pair of M g {i) and y max , and extract from it our theoretical prediction 
and uncertainty Ai?^ h for these parameters. We then take the ratio R^/R^, combine the 
theoretical and experimental errors in quadrature, and test at which confidence level this ratio 
is compatible with 1. In both versions of the RS model, the data exclude significant portions 
of the model parameter space. With the conventional choice ?/ max = 3, for example, one finds 
M (i) > 13 TeV for the brane-Higgs model and M.(i) > 4.5 TeV for the bulk-Higgs scenario, 
both at 95% CL. Weaker constraints are obtained for smaller values of y max . These bounds 
may be compared with those derived from the analysis of electroweak precision observables. 
The strongest constraint arises from the S and T parameters [56], whose present values are 
S = 0.03 ± 0.10 and T = 0.05 ± 0.12, with a correlation coefficient p = 0.89 [57]. In the 
minimal RS model, one obtains at tree level [58] 

2irv 2 ( 1 \ „ ttv 2 ( „ 1 



S Ml K V 1 L ) ' T 2 cos 2 9 W Ml K \ L 2LJ- ( " } 

Requiring that these corrections are compatible with the experimental data, we find that 
M g (i) > 12 TeV at 95% CL. This strong bound, which is indicated by the dashed line in Fig- 
ure 4, may however be weakened in several ways, for instance by including loop corrections, by 
reducing the size L of the extra dimension (so-called "little RS models") [59], or by introducing 
large brane-localized kinetic terms in the RS Lagrangian [58]. We note that for M (i) > 12 TeV 
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M gW [TeV] M gW [TeV] 

Figure 5: Predictions for the ratio Rh in the custodial RS model, for the cases of a 
brane-localized Higgs boson (left) and a narrow bulk-Higgs field (right). The meaning 
of the colors and curves is the same as in Figure 3. 



there is no significant flavor problem of the minimal RS model, as the tightest constraint from 
the ex parameter in K-K mixing [9] can be satisfied with a modest 25% fine-tuning [60]. 

Softening the constraints from electroweak precision tests by means of a symmetry is the 
main motivation for extending the RS model by enlarging the gauge group in the bulk [26- 
28]. We now proceed to study the RS model with custodial symmetry, in which the Wilson 
coefficients C\ and C5 are given in (94). The corresponding numerical results are shown in 
Figure 5. For large masses M g (i) we can derive analogously to (98) a formula for Rh depending 
explicitly on N g and 2/ max , which in the present case reads 



Rh 



v 2 



2M 2 



V 2 



(100) 



Note that the leading terms proportional to N 2 are enhanced by a factor 4 compared with the 
minimal model, reflecting the larger multiplicity of KK quark states. The remaining terms 
are enhanced by a factor 2, as can be seen from (92) and (96). As a result, in the custodial 
RS model one finds significantly larger corrections to the SM prediction Rh = 1 than in the 
minimal model [22]. The terms in brackets then evaluate to approximately [78.7 y"^^ + 7.1] 
for the RS model with a brane-localized Higgs, and [— 65.3£/^ ax + 7.1] for the model with a 
narrow bulk Higgs. 

Figure 5 confirms the fact that the corrections to the Higgs-boson production rate are 
much enhanced compared with the case of the minimal RS model. Correspondingly, we obtain 
significantly larger exclusion regions than for the minimal model. This is shown in Figure 6. In 
the brane-Higgs scenario, we obtain the exclusion range 4.5 TeV < M g {v, < 19 TeV for y max = 3 
at 99% CL, while in the bulk-Higgs model we find the lower bound M g m > 9.5 TeV at 95% CL. 
Note that the allowed region in the upper left corner (at small M g m and large y max ) of the first 
plot in the figure is one in which the new-physics contribution to the gluon fusion amplitude 
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Figure 6: Excluded regions of parameter space in the custodial RS model, for the 
cases of a brane-localized Higgs boson (left) and a narrow bulk-Higgs field (right). The 
vertical dashed line shows the lower bound on M g {v, obtained from a tree-level analysis 
of electroweak precision observables. 



is larger than the SM contribution by about a factor 2 and interferes destructively, which 
appears somewhat unnatural. Moreover, it has been argued that most models in which the 
gluon fusion amplitude has the opposite sign than in the SM have problems with fine-tuning 
and vacuum stability [61]. The bounds on the RS parameter space that can be derived from 
Figure 5 are stronger than those derived from the analysis of electroweak precision observables. 
In the custodial model the formula for the S parameter shown in (99) is left unchanged, while 
the custodial protection removes the leading term proportional to L in the expression for the 
T parameter, such that T = —irv 2 / (4L cos 2 Ow^kk) P6]- Requiring that these corrections are 
compatible with the experimental data, we find that M (i) > 4. 7 TeV at 95% CL. As indicated 
by the dashed line in Figure 6, this lower bound is generally much weaker than the constraints 
implied by Higgs physics, except for regions in parameter space where y max is very small. 
Note that for such small values of the KK mass scale but ?/ max ~ 3, the RS flavor problem for 
the ex parameter can be solved by a fine-tuning of 5-10%, or alternatively by enlarging the 
strong- interaction gauge group in the bulk [60]. 

We may also read the exclusion plots in a different way. If we would like to have the first 
KK excitations in the reach for direct production at the LHC, then this imposes a strong 
upper bound on the maximum allowed values of the elements of the 5D Yukawa matrices. For 
instance, assuming that M g w = 5 TeV, we find that y ma x < 0.6 in the brane-Higgs model, and 
2/max < 1-5 in the bulk-Higgs scenario (both at 95% CL). Too small Yukawa couplings would 
however give rise to enhanced corrections to €k [9], and hence they would reinforce the RS 
flavor problem. 

The above analysis shows that Higgs physics, and in particular the Higgs-boson production 
rate in gluon fusion, provide sensitive probes of the virtual effects of KK excitations in the 
context of various RS scenarios. While models with a brane-localized scalar sector predict a 
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M gW [TeV] M gW [TeV] 

Figure 7: Predictions for the ratio Rh in type-II brane-Higgs models with two different 
Yukawa matrices in the couplings of the Higgs doublet to Z2-even and Z2-odd fermion 
fields, as shown in (C.l), for the cases of a brane-localized Higgs boson (left) and a 
narrow bulk-Higgs field (right). The meaning of the colors and curves is the same as 
in Figure 3. 



suppression of the gluon fusion rate, this rate tends to be enhanced in scenarios with a bulk- 
Higgs field. The two classes of models can thus easily be distinguished in their signatures. 
The bounds on the model parameters obtained from Higgs physics are complementary to and 
sometimes stronger than those derived from the analysis of electroweak precision observables 
and rare flavor- changing processes. In models with a custodial protection of electroweak pre- 
cision observables, the indirect effects of KK states on the Higgs-boson production rate are 
strongly enhanced compared with minimal RS models, and hence Higgs physics provides the 
strongest constraints in this case. Even under the pessimistic (but not unrealistic) assumption 
that the direct detection of KK excitations is out of the reach of the LHC, one may still see 
sizable modifications of the Higgs-boson production cross section. For example, even with 
M g (i) = 10 TeV or even 15 TeV, Figures 3 and 5 show that virtual effects of KK particles 
can have significant effects on the Higgs-boson production cross section, provided that the 5D 
Yukawa couplings are not too small. We also note that different implementations of warped 
extra-dimension models, such as little RS models in which the bounds from electroweak preci- 
sion measurements and flavor physics are relaxed by reducing the size L of the extra dimension 
[59], give rise to very similar Higgs phenomenology, because the L-dependent corrections in 
(98) and (100) have only a minor impact. 

We close this section with a brief investigation of the variant of the brane-localized Higgs 
scenario with two different, random Yukawa matrices and Y q s ^ in (8), which was discussed 
in Section 5.7 (see also Appendix C). We treat these matrices as completely independent and 
impose the same upper bound ?/ max on their entries. We refer to these scenarios as type-II 
brane-Higgs models. With the exception of the effect of the shift of the Higgs vev encoded by 
k v , the new-physics effects arising in this model have indefinite sign. Also, the leading terms are 
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no longer enhanced by a factor Ng = 9 as in (98) and (100), but rather the sum of Ng random 
complex numbers of variable phase scales like N g (random walk). The resulting corrections on 
the production rate are therefore smaller in magnitude. This is shown in Figure 7, which shows 
our results for in both the minimal (left) and the custodial RS model (right). Because 
of the distribution of points around a central value close to the SM prediction = 1, it 
is impossible to exclude any regions of parameter space in type-II brane-Higgs models. On 
the other hand, one may still find appreciable effects on the Higgs-boson production rate in 
these models. For example, with y max = 3 a modification of Rh by 20% is possible for KK 
excitations as heavy as lOTeV in the minimal RS model, and 13TeV in the custodial model. 

8 Conclusions 

The discovery of a Higgs-like boson at the LHC in the summer of 2012 [1, 2] has marked the 
beginning of a new era in particle physics. The demand for an explanation of the hierarchy 
problem is now more pressing than ever. Precise measurements of the Higgs-boson couplings to 
various SM particles can provide valuable tools to distinguish between different new-physics 
models, which can address the hierarchy problem. Such measurements can elucidate the 
mechanism of electroweak symmetry breaking and probe for indirect hints of new particles. 
Of particular importance are loop-induced processes, such as the Higgs-boson production 
process gg — >■ h and the radiative decay h — > 77, since possible new heavy resonances can lead 
to sizable deviations from the SM expectations. 

In this paper, we have focused on the gluon fusion cross section in various incarnations 
of RS models, in which the scalar sector is localized at or near the boundary of a warped 
extra dimension. We have derived an exact expression for the gg — )■ h amplitude in terms of 
an integral of the mixed-chirality components of the 5D fermion propagator with the Higgs- 
boson profile along the extra dimension, given in (15) -(17). This expression can be used 
to calculate the effective CP-even and CP-odd ggh couplings defined in (9), as long as one 
succeeds in deriving an explicit expression for the propagator. In contrast to the procedure 
commonly used in the literature, all our calculations have been performed by keeping the exact 
dependence on the Higgs-boson mass. Moreover, working in a 5D framework, we have not 
distinguished between the SM particles (zero modes) and their KK excitations. In our final 
results for the effective ggh couplings, however, we have been able to identify the contributions 
of the SM particles and the infinite tower of heavy KK states. 

In Table 1, we have classified different versions of RS models according to the parametric 
relation of the characteristic width rj of the Higgs-boson profile with respect to the two ratios 
v\Y q \/M-^yi and v\Y q \/AiceY, where AT e v is the value of the inherent UV cutoff near the IR 
brane. We have shown that it is possible to obtain explicit analytic expressions for the 5D 
propagator for both a brane-localized scalar sector and a scalar sector that lives very near 
the IR brane (narrow bulk-Higgs scenario). To an excellent approximation, the effective ggh 
couplings in these cases only depend on the 5D Yukawa matrices Y q and the ratio v 2 /M^ K , see 
e.g. (70). On the contrary, the results for a generic bulk-Higgs model, in which the width of 
the Higgs profile is not parametrically small, depend in a complicated way on the 5D fermion 
masses and the shape of the Higgs profile [44] . 
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Importantly, we have pointed out that there is no controllable interpolation between bulk- 
Higgs and brane-Higgs models. In RS models in which the scalar sector is localized on the IR 
brane, one finds that the gluon fusion cross section is reduced compared with its SM value, 
in accordance with the findings of [19, 22, 23]. In this context, we have proved a conjecture 
made in [23] for the analytic form of the contribution from virtual KK states. On the other 
hand, in models in which the Higgs-boson is described in terms of a narrow bulk field localized 
near the IR brane, the cross section is enhanced (apart from regions in parameter space in 
which the 5D Yukawa matrices have very small entries). This result confirms the calculations 
performed in [20]. The qualitative difference between the predictions obtained in the two 
types of scenarios provides an opportunity to distinguish between the two classes of models, 
provided that a deviation of the gluon fusion rate from its SM value is observed in the future. 
When one tries to interpolate between the bulk-Higgs and brane-Higgs scenarios, for instance 
by considering the limit rj — > in the context of a bulk-Higgs model, one enters a transition 
region with 77 ~ v\Y q \/A^ e Y, in which the contributions from higher-dimensional operators in 
the RS Lagrangian become unsuppressed or even enhanced, so that the effective field-theory 
approach breaks down. 

For most of our discussion, we have focused on the minimal RS model with gauge symmetry 
SU (3)c x SU (2)^ x U (l)y in the bulk. However, in Section 6 we have considered an extension 
with a custodial symmetry, based on the gauge group SU(3)c x SU(2)l x SU(2)r x U(1)x x 
Plr- We have succeeded in deriving analytical expressions for the effective ggh couplings in 
terms of the same input parameters that appear in the minimal model, see (94). We have found 
that the contribution from the infinite KK tower of virtual quark states is four times larger 
than in the minimal model, due to the higher multiplicity of particles running in the loop. We 
have also investigated a variant of the brane-Higgs scenario (dubbed type-II models), in which 
one introduces two different 5D Yukawa matrices coupling the scalar doublet to Z 2 -even and 
Z 2 -odd fermion fields. 

Investigating the phenomenological implications of our results, we have focused on the 
ratio Rh representing the gluon fusion cross section in the various RS models normalized to 
its SM value. We have distinguished between a brane-localized and a narrow bulk-Higgs field 
for both the minimal and the custodial RS model. We have pointed out the fact that the KK 
contribution to Rh does not only depend strongly on the number of quark generations (N g = 3), 
but also on the maximal value y max one imposes on the magnitudes of the individual entries 
of the anarchic 5D Yukawa matrices, which are assumed to be random complex numbers. 
To a good approximation, our results can be parameterized in terms of ?/ max and the mass 
M g (i) ps 2.45Mkk of the lightest KK gluon state. Provided that the value of y max is not too 
small, we have shown that quite generically the new-physics effects in RS models can lead 
to significant deviations of Rh from 1, even for KK masses that are not in the reach of the 
LHC. This is evident from Figures 3 and 5, which show that significant corrections can be 
obtained even for KK gluon masses in the range of 10-20 TeV. For RS models with a custodial 
symmetry, whose original motivation was to lower the KK scales via a protection of the T 
parameter and the Zbb vertex, the effects are even more pronounced. As mentioned earlier, Rh 
is strictly less than 1 in RS models with a brane-localized scalar sector, whereas it exceeds 1 
in models with a bulk-Higgs field for almost all points in parameter space. 

Comparing our predictions with the latest ATLAS and CMS data [55], we have derived 
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exclusion regions in the M g (i) -y max parameter space of the various models. The corresponding 
results shown in Figures 4 and 6 demonstrate the new-physics reach of Higgs-boson observables 
such as Rh in an impressive way. In the minimal RS model and at 95% CL, one can exclude KK 
gluon masses lighter than 12.8 TeVx (y max /3) for the brane-Higgs case and 4.4 TeVx (y max /3) 
for the case of a narrow bulk-Higgs field. In custodially protected RS models, these bounds 
increase to 24.4 TeVx (y max /3) and 9.6 TeVx (?/ max /3), respectively. Especially in this latter 
case, the bounds derived from Higgs physics are already much stronger than those obtained 
from electroweak precision tests. A possible way to weaken these bounds is to assume that 
2/max is significantly smaller than the value commonly adopted in the literature (y max ~ 3). 
However, this would create a tension with other observables, such as the parameter €k in the 
neutral-kaon system, which in the context of RS models receives corrections scaling like 1 / y^ ax 
[9]. Another possibility is to focus on type-II brane-Higgs models, see Figure 7. We have shown 
that in these models the new-physics effects on the gluon fusion cross section are not as large 
as in RS models with only one Yukawa matrix Y q , and due to an undetermined complex 
phase they can interfere both constructively and destructively with the SM contributions. For 
l/max ~ 3, this can lead to an enhancement or a suppression of the cross section, which can be 
in the 20% range even for KK excitations as heavy as 10 TeV in the minimal and 13 TeV in 
the custodial RS model. 

The methods developed in this work can be extended to analyze the loop-mediated decay 
h 77 as well as other decay modes of the Higgs boson. As the experimental precision on the 
extracted Higgs couplings increases, it will be exciting to confront the theoretical predictions 
obtained in various RS models with the data. One might hope that, perhaps, one day this 
could provide a first hint of the possible existence of a warped extra dimension, even if no KK 
excitations of SM particles are to be discovered at the LHC. 

Note added: While this paper was in writing, the work [62] appeared, in which similar 
questions as in the present work were addressed. While we have no objections to the ana- 
lytical calculations presented in that paper, we disagree with the interpretation of the results 
obtained by these authors. In particular, the argument that higher-derivative operators in 
the RS Lagrangian would "dress" the brane-localized Higgs to make it look like a bulk field 
is incompatible with our findings. Rather, these operators dress the bulk Higgs as its profile 
is made narrower, and they are responsible from the transition from an enhanced gg — >■ h 
amplitude (bulk Higgs) to a suppressed one (brane-localized Higgs), see Table 1. 
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A Details of the solution for the propagator functions 



Here we present in more detail the derivations needed to calculate the functions T±(p E ) defined 
in (15). Since ultimately we only need the mixed- chirality components of the 5D fermion 
propagator evaluated with t = t' and convoluted with the profile of the Higgs boson along the 
extra dimension, we can from the beginning assume that 1 — 77 < t' < 1, but we allow t to 
take any value. 



Calculation of the propagator functions A. LL and A RL 

For t < 1 — 77 in the bulk, the most general solutions for the propagator functions are superpo- 
sitions of modified Bessel functions, as shown in (31) for the case of A LL . The function A RL 
then follows from the second equation in (23). Imposing the boundary conditions (25) on the 
UV brane (at t = e), one finds four relations among the right coefficients Ci(t'). Rescaling 
these coefficients appropriately, we write the solutions in the form 



,., 2> r t (D?(p E ,t) \ I K,(t) K 2 (t') 

where (with A — Q,q) 

Di(pE,t) = I_ CA i(ep E )I CA i(pEt) - I CA+ i(ep E ) I_ c . + i(p E t) 



(a.i; 



D£(p E ,t) = I_ CA _i(ep E ) I CA+ i(p E t) - I CA + i{ep E ) I_ c i(p E t) 



(A.2) 



are diagonal matrices, and £)^(p£,e) = 0. 

In the region t > 1 — 77 very near the IR brane, the general solution for A q LL (t, t'; —p 2 ) has 
been given in (32), while the second equation in (23) yields 



A * ( t t > _ 2x ^kk / S q S(t) gY q C(t) \ / d(f) C 2 (t') 
' V \QX}C(t) S q S(t) J \C 3 (t>) C A (t>) 

+ M KK ( S q C{t) gY q S{t)\ ( C 5 (t>) C 6 (f) 



(A.3) 



V \gYjS(t) S q C(t) J \C 7 (t>) C 8 (t>) 

where g = u/(v2Mkk)) and we have used the abbreviations 

C(t) = cosh[5 g 0"(t - 1)] , S(t) = sinh[5 g 6 r '(t - 1)] , (A.4) 

and similarly C(t) and <S(t) defined with S q instead of S q . Because of the discontinuity at 
t — t', we must distinguish the cases where t > t' and t < t'. We indicate this by means of a 
superscript on the coefficient functions, using the notation Cf{t') for t > t', and Cf(t') for 
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t < t' . Imposing the boundary conditions (25) on the IR brane (at t — 1), and satisfying the 
jump conditions (24), it is straightforward to show that 

c>(t') = c<{t') + x3— ^-5(0 , c<(t') = A-^-m , 

KK °9 KK J 9 

C 2 >(f) = £<(*') , C>(t') = <%<(*') , (A.5) 

JW KK JW KK 



while all other coefficients vanish. These relations allow us to express the solution in terms of 
the four functions Cf(t') with i = 1,2, 7, 8. 

The remaining eight coefficients are determined by requiring that the solutions for the 
propagator functions be continuous at t = 1 — rj. Continuity of A LL yields the conditions 



D?(Pe,1 


-ri)K x {t') 


= cosh S q C<{t') + 




- V )K 2 (t') 


= cosh S q C<(t'), 






= smhS q Cf(t') , 




- V) K 4 {t') 


= smkS q C<(t') - 



(A.6) 



77 COsh Sq 



S(t>) 



which can be used to eliminate the coefficients Ki(t'). Note that on the left-hand sides of these 
equations we can take the limit rj — > without difficulty. When the solutions are inserted into 
the expression for Are in (A.l), we then encounter the ratios Ra(Pe) = D^(pE, 1)/ D^ipE, 1) 
defined in (35). The remaining four coefficients Cf(t ; ) can be derived by requiring that the 
propagator function A_r£ is continuous at t — 1 — r\. To express the answers in a compact 
form, we introduce the definitions 



Nf(p%) = 1 + Zf{pl) + v p E [l + R Q \p E ) (l?)- 1 R q (P E ) I?] R Q (Pe) 



Nf{p 2 E ) = 1 + Zf{p E ) + V p E [R q (Pe) + Y q R q (p E ) Y- 1 ] 



(A.7) 



where 



S 9 



KK A g J 9 



(A.8) 



In the limit 77—7-0, the quantities N^'*(p E ) approach 1 + Z q {p 2 E ) with Z q {p 2 E ) as defined in 
(52), while the quantities Z^' % {p E ) approach Z q (p 2 E ). After some lengthy algebra, we now 
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obtain 



cm 
cm 
cm 



M£ K S q smhS q 



x 



sinh 2 S q + Zf(p 2 E 



Nf(pl) 



Rqipi 



We 



1 



1 



PeM kk cosh S q 

1 



PeM kk S„ cosh 5, 



Rq(Pe 



1 5(f) 
^ ^* 2 (P|) S q coshS q 



9 • 



Nf{pl) cosh 5,' 



M| K XfsinhS, 



x 



sinh S„ 



Nf{p 2 E ) 



cosh S„ 



Nf(pl) 



SgtanhSg 

— — : RqVPl 

VPe 



Nf{p 2 E ) 



Sit') 



S q COSh Sq 



gYq. 



(A.9) 



Calculation of the functions T±(p E ) 

Equipped with all required coefficients, we can now derive explicit expressions for the quantities 
T±(p E ) defined in (15). Using (A. 3), we find that 

JL( Y q \ A^ L (M;p|) + AUM;p|) 
V2\YJ J 2 



2rj 



C{t)C<{t) 



S{t)C<{t) 



V S(t)C{t) 



M 2 KK 



+ gYjS g S(t)C<(t) 



+ QY q SqC(t)C<(t) +h.c.\, 



(A.10) 



where the contribution from A Q LR is the hermitian conjugate of that from from A^, L (assuming 
p\ > for now). Upon taking the trace in (15), the two terms proportional to rj in the square 
brackets cancel each other. Next, using the explicit expressions for the coefficients in (A.9), 
we find that the contribution involving the terms proportional to I/77 in square brackets in 
the expression for Cf cancel against the contribution from in (A. 10), and likewise for the 
terms involving and Cf. After the dust settles, we obtain 



T + ipl) 



£ fdtSl{t-l)T, 

q=u,a 



S q sinh 2S q 



x 



sinh 2 S q + C 2 {t)Zf{p% 



N^ipi 



~ * W » 7-77,2/ 2 \ 



(A.11; 



+ h.c. 
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and analogously for T_(p 2 E ). In the case of one generation, the above expression reduces to 
formula (42), once we identify k\{p E ) = 1 + Z q {p 2 E ) coth 2 S q . For the general case of three 
generations, the result (A. 11) simplifies if we take the limit rj — > 0, in which we recover the 
results shown in (51). Note that in this case the dependence on t inside the square brackets in 
(A. 11) disappears, due to the identity C 2 (t) — S 2 (t) = 1. Therefore, as already mentioned in 
Section 5.5, we would have obtained the same result by setting t = t' = 1~, as shown in (56). 



Generalizations for the model with custodial symmetry 

The derivation of the propagator functions in the RS model with custodial symmetry proceeds 
in an analogous way. In fact, the only difference arises in the equations in (A.l), where Di 2 
and D\ 1 must be replaced by 



D? t2 (p E ,t) 

D®(p E ,t) 



(Df,{p El t) 



and 



V 



o 

Dl] 3 (p El t) 
Dl%(p E ,t)J 



(A.12) 



for up-type quarks, and analogously for down- and A-type quarks, with patterns that can be 
read off from (85). The appearance of the functions 



-CA- 



-CA- 



.{ep E ) 1 



CA 



ipEt) 
AVEt) 



CA~ 



CA~ 



{ep E ) I_ CA+ i(p E t) 



(A.13) 



gives rise to the ratios R A \p E ) = D^(p E , l)/Df(p E , 1) defined in (82) 



B Case of a bulk-Higgs field 

This section intents to relate an RS model with a scalar sector in the bulk, in which the Higgs 
field and its vev have profiles that are strongly peaked near the IR brane, to the RS model 
with a brane-localized Higgs sector. Our discussion will follow the expositions given in [38, 63], 
but we will generalize these results in some aspects. 



Definition of the model 

Using the orbifold coordinate 25 = r<f), the action for the Higgs sector reads 

d 4 x / dx 5 e- 4 ^ g MN D M ^D N ^-fi 2 \^\ 2 -V VY ^)S(x 5 )-V m (^)S(\x 5 \-r7r) 

J —rir 

(B.l) 

where \x provides a bulk mass for the scalar field, which can be tachyonic (see below). The 
potentials localized on the UV and IR branes determine the boundary conditions of the scalar 
fields and induce electroweak symmetry breaking. They are chosen to be 

Vuv($) = M uv |$| 2 , V m ($) = -M m |$| 2 + Air |$| 4 , (B.2) 
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with mass dimensions [Muv] = [Mm] = 1 and [Air] = —2. The dimensionful parameters in the 
5D action naturally scale with appropriate powers of Mpi, and we find it useful to introduce 
dimensionless 0(1) parameters by the rescalings 

M uv M m X m k 

muv = ^T' miR "^F' A "^T- (R3) 

We now change variables from to t — e e ff ^' and express the scalar doublet in 
the form 

e v^ V75 [«(*) + t) + ifo(x, t)])' 1 ' ; 

where v(t) denotes the profile of the Higgs vev along the extra dimension, h(x,t) is the 
5D physical Higgs scalar after electroweak symmetry breaking, and <f) + (x,t), 4>s(x,t) are 5D 
Goldstone bosons. For the following analysis we do not consider the Goldstone fields any 
further (unitary gauge). Integrating by parts, the Lagrangian corresponding to the action 
Sh = J d x Ch(x) in (B.l) can be rewritten in the form 



£h(x) 



+ 



Mkk 



V(t) + 2H(X, t) + _ fl tg + l|M {t 2 d 2 + Wt _ 02) KX, t) 



t V 1 ' t t V 1 ' ' t 

-i r 

d t [tv{t)] + '-^f^d t [th{x,t)} 



7rM| K rft;(t) + 2/i(a;,t) h(x,t) 



p 2 



[w(e) + /i(x,e)] : 



- mm [v(l) + /i(x, l)] 2 + [v(l) + l)] 4 



(B.5) 

where /3 = \J & + /i 2 //c 2 . Requiring that the terms linear or quadratic in /i(x, £) cancel on the 
UV and IR branes yields the boundary conditions 11 

2A 

dt[tv(t)} t=e+ = m vv v{e), d t [tv{t)] t=1 . = m m v(l) - — — v 3 (l) , 

[t /i(ar, t)] t=e+ = m uv e) , d t [t h(x, t)] t=1 _ = m m h(x, 1) - — r - v 2 (l) h(x, 1) . 

KK 

The notation e + and 1~ indicates that the orbifold fixed points must be approached from the 
appropriate sides. 

Profile of the Higgs vacuum expectation value 

By means of the variational principle with respect to v(t), one obtains the equation 

(t^ + td t -(3 2 )^- = 0, with /3 2 =4 + ^ (B.7) 
t k 



11 These conditions can also be derived by integrating the field equations over infinitesimal intervals about 
the branes. 
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which ensures that the tadpole terms in the Lagrangian (B.5) cancel out. We then obtain 



Ch(x) 



2tt 
~L 



1 dt 

T 



7T 



A 



\ d,h(x, t) d»h(x, t) + ^ *M (t 2 d 2 + td t - /3 2 ) *M 



v 4 (l) + Av(l)h 3 (x,l) + h 4 (x,l) 



(B.8) 



The general solution to the differential equation (B.7) subject to the boundary conditions 
(B.6) is 

2 2 + - m uv 



and 



i(t) = N v - r v t l ~ p ) , with r. 
(mm - 2 



N, 



2A 



2 - f3 - m uv 
P) ~ r v (m IR -2 + (3) 



1 



(B.9) 



(B.10) 



Before proceeding, let us first discuss which values the parameter /3 can take. Motivated 
by the observation that the energy-momentum flux in a pure anti-de Sitter space without an 
IR brane (which corresponds to taking the limit r — > oo) vanishes at the boundary only if the 
5D scalar field obeys the Breitenlohner- Friedman bound fi 2 > —4k 2 [64], one usually assumes 
that (3 must be a real positive number, even though not necessarily larger than 2. Unless (3 is 
very close to zero, it follows that the coefficient r v oc e 2/3 in (B.9) is extremely small and can 
be set to zero for all practical purposes. The only exception would be the region where t ~ e 
is very near the UV brane, which however is irrelevant for our analysis here. It follows that 



v(t) =v(l)t 1+l3 



with v(l) 



Mi 



KK 



m m -2-/3 



2A 



(B.11) 



The requirement that the Higgs vev be a real number imposes an upper bound on the param- 
eter (3, since A > is required by vacuum stability. We thus obtain the allowed range 



< (3 < m m - 2 



(B.12) 



We proceed to relate the parameter v(l) to the physical value vsm of the Higgs vev in 
the SM. After electroweak symmetry breaking, the mass terms for the W and Z bosons are 
generated by the 5D Lagrangian 



d A x 



2tt [Ut v 2 (t)g 2 



W+{x,t) W-i*(x,t) + 



1 



2 cos 2 6 



Z^x,t) Z"(x,t) 



(B.13) 



where the 5D gauge coupling g 5 is related to the gauge coupling g of the SM by g = g^jyhvr- 
\. Introducing the KK decomposition 



^ oo 



(B.14) 



n=0 
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and similarly for the W bosons, and using that the zero-mode profiles are flat, Xn W = l/v^vr 
up to higher-order terms in v 2 /M^ K [65], we can identify 



_ 27T f 1 dt 2 7T v 2 (l) 



* S T). T^ = ITTi- (B15) 

where once again we neglect terms suppressed by powers of e. It follows that 

v(t) = v 4 + . (B.16) 

The parameter V4 coincides with the parameter v used elsewhere in this paper. At lowest order 
in an expansion in powers of v 2 /M^ K , it coincides with the SM parameter vsm as defined, e.g., 
via the value of the Fermi constant. Higher-order corrections to the relation v$m = v 4 could be 
calculated by solving the differential equations for the profiles of the gauge-boson zero modes 
in the presence of the Higgs vev. 

Profiles for the Higgs boson and its KK excitations 

We now proceed to study the eigenvalue problem for the physical Higgs boson and its KK 
excitations. We write the KK decomposition of the 5D Higgs field as 

oo 

h(x,t)=^T l h n (x)xn{t), (B.17) 

n=0 

where the zero mode h(x) = ho(x) corresponds to the SM Higgs boson. The profile functions 
obey the orthonormality condition 

2tt [ x dt 



^ J j Xm(t) Xn(t) = 5 mn , (B.18) 



which ensures that the kinetic terms in the effective 4D Lagrangian are canonically normalized. 
In order to obtain canonical mass terms from the Lagrangian (B.8), we must impose the 
equation of motion 

(t 2 d 2 t + td t + t 2 x 2 n - (3 2 ) ^ = , (B.19) 

where x n = m n /MKK denote the masses of the KK scalar bosons in units of Mkk- The general 
solution to this equation is a linear combination of Bessel functions, 

X n(t) = N n t [Mx n t) - r n Yp(x n t)] , (B.20) 

where the boundary condition on the UV brane in (B.6) once again implies that r n oc e 2/3 is 
extremely small and can be set to zero for all practical purposes, since we are not interested 
in the region where t ~ e. We then obtain 

Xn(t) = , tMXnt) • (B.21) 
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The boundary condition on the IR brane gives rise to the eigenvalue equation, which deter- 
mines the masses of the scalar modes. We find 



2(m m -2-/3) =26. 



(B.22) 



It follows from this equation that even the zero mode (the SM Higgs boson) has a mass that is 
naturally of order the KK scale Mkk, which empirically cannot be less than a few TeV. This 
is the little hierarchy problem, which as mentioned in the Introduction is not addressed in RS 
scenarios. In order to obtain a realistic Higgs mass <C Mkk, we must assume that 



5 = m m - 2 - < 1 . 



(B.23) 



Once this is done, it is straightforward to obtain a formula for the zero-mode mass in a power 
series in 5. We find 



mi 



■I'r 



M 2 KK 



4(1 + 0)5 



+ 



25 2 



+ 



(B.24) 



2 + P ' (2 + P) 2 (3 + 0) 

Assuming Mkk = 2 TeV, for example, implies that (1 + P) 8 ~ 10~ 3 , which corresponds to 
a fine-tuning of 1 in 1000. For the zero-mode profile, it is now straightforward to obtain an 
expansion in powers of Xq. The leading terms are given by 



Xo{t) 



L 



71 



(1 + 0) 



X n 



t 2 



p 



p 



(B.25) 



Dropping the irrelevant constant proportional to i> 4 (l), the Higgs Lagrangian (B.8) can 
now be written as 



-d^h n (x) d^h n (x) 



m 



n 1.2 



h n (x) 



4L 



-v A 



L 



7Y 



1 + Pf\J2 6 £rn £n h{x) h rn (x) h r 



[X 



(B.26) 



l,rn,n 



— (1 + /3) 2 A y~] £m€nh k (x)hi(x)h m (x)h n \ 

71 , +^ 



X 



k,Lm,n 



where £ n = Xn(l)/y ~ (1 + P)- From (B.25) we find £ ~ 1 f° r the zero mode, while (B.21) 

and (B.22) imply that ~ ±1 / a/I + ~j3 for the KK excitations. We proceed to relate the 
parameter A to the physical value A4 of the Higgs self coupling. The relevant terms in the SM 
Lagrangian are 



in 



h 7,2 



^sm Asm h 



Asm h , 



(B.27) 



2 u ° 4 

where m\ = 2Asm^sm- Matching either one of these terms with the corresponding term in 
(B.26), we obtain at leading order 



4L 

Asm = A 4 = — (1 + P) 2 A = A IR k 2 (1 + P) 2 . 

71 



(B.28) 
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The relation between Asm and A4 receives higher-order corrections in v 2 /M^ K , which are 
calculable in the model and depend on which of the three couplings in (B.27) is used to 
perform the matching. 



Yukawa interactions 

We finally consider the Yukawa couplings of the scalar field to the fermions. In the model 
with a brane-localized Higgs sector, which we have considered for most of this work, one has 
in analogy with (6) 



^brane | 



X 



q=u,d 



V'^ f - 1 > + ™- 1 W)- 



V2 



5D 



Y q 




Q R {t, x) + h.c. , 
(B.29) 

where the 5D Yukawa matrices Y^ D have mass dimension —1. In the model with a bulk-Higgs 
field, we have instead 



-L 



bulk J 



q=u,d 



1 dt v(t) + E n K( X )x n (t) e 1 



y/2 





1 <?,bulk 



-5D \ 

9 Q bulk J Q R (t,x) + h.C. 



(B.30) 

where the 5D Yukawa matrices Y^^ now have mass dimension —1/2. In order to match 
the two expression onto each other, we must rewrite the functions v(t) from (B.16) and Xo(t) 
from (B.25) in terms of functions with unit area, which can be mapped onto the normalized 
distributions 5%(t — 1) and <5^(t — 1). We obtain 



with 



71 2 + f3 



Xo(t) 



71 2 + P 



1 + 



4(l + /3)(2 + /3)(4 + /9) 



+ 



81^^-1) = (2 + P)t 1+ P 



8 1 h / \t-l) = (2 + P)t 



5}!\t - 1; 



p 



(B.31) 



(B.32) 



4(1 + 0) 

Here 1/(3 plays the role of the regulator 77 in (B.29). Using the quark bilinear terms as a 
reference, the corresponding matching relations between the two Yukawa matrices reads 



k 



Y — — Y 



5D 



Vk(l + (3) 
2 + /3 



v 5D 

^g.bulk 



(B.33) 



The quantities on the left-hand side of the equation are the dimensionless Yukawa matrices 
introduced in (6), whose elements are assumed to be random numbers bounded in magnitude 
by 2/ m ax- If one used the hqq couplings instead, the above relation would receive corrections of 
0{xl). 
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Limit of a narrow bulk-Higgs field 



We are finally in a position to study the limit /3 3> 1, in which the profile functions in (B.32) 
become strongly localized near the IR brane. We can then identify 1/(3 with the width of 
the Higgs profile, which plays the role of the regulator 77 in our brane-Higgs scenario. The 
Yukawa matrices of the bulk-Higgs model must then be identified with Y q -h- \/ k / /3 Y q 5 ® ulk = 
(k/y/Ji) Y^buik- Finally, since t is pushed near 1, we conclude from (B.32) that 

as was claimed near the beginning of Section 4. 

Taking the limit of very large f3 is not particularly natural, since IS 
naturally of 0(1). For large /3, we have the double hierarchy 

— <C k <C u ~ — — , or — ~ k ~ — . (B.35) 

r 2 r p 

Large /3 can be achieved by taking k significantly smaller than the Planck scale (and 1/r yet 
smaller by an order of magnitude), or by assuming that /1 and Mir are significantly larger 
than Mpi. The first possibility appears more plausible. Note that for large /3 relation (B.28) 
implies that A4 ~ Air/i 2 , indicating that increasing /3 by lowering the curvature parameter k 
does not affect the relation between A 4 and Air in a significant way. 



C Case of two different Yukawa matrices 

We briefly discuss the generalization of our results to the case where the two Yukawa couplings 
in (8), involving products of Z2-even and Z2-odd fermion profiles, are associated with different 
Yukawa matrices, such that 

C n = ^ fdtnt - 1) [ogOt C £\t) Y? CM® a?' + a£>t S«>(t) Y u s ^(t) a$0 

(C.l) 

At the level of the gluon fusion amplitude (12), the above modification is implemented by the 
substitution 

iv Y » ) " ^ * y* + ^ (i?* y o ) pl ■ <a2) 

This generalization is only allowed if the Higgs boson is localized on the IR brane. For a 
bulk-Higgs field, it is forbidden by 5D Lorentz invariance, since 275 is one of the 5D Dirac 
matrices 7°. 

The equations of motion (23) for the propagator functions must now be generalized to 



p 2 AL(M'; -P 2 ) - M KK (I + M q (t)) AUtA -P 2 ) = 5{t - 

a^(m';-p 2 )-Mk K (-| + -^5w) K L (t,t';- P 2 ) = o, 



(C.3) 
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where 



M q (t) = U 



C Q 




V2M K 



w - 1; 



KK 







yc 




(C.4) 



replaces the generalized mass matrix in (20). The coupled set of first-order differential equa- 
tions in (C.3) can be combined to yield the second-order equation 



^-M q {t)M\{t) 



dM\(i) 
dt 



Pe 



1 



-p 



M 2 KK 



8{t-t') 



(C.5) 

In the bulk region t < 1 — 77, where the profile 8% (t — 1) of the Higgs vev vanishes and the 
mass matrix is hermitian, this equation reduces to the original equation (26). However, its 
structure becomes much more complicated for t > 1 — 77. We have not succeeded to derive the 
general solution in that region. 

In the case of infinitesimal 77 (at fixed p 2 ), however, it is consistent to only keep the terms 
in (C.3) that are enhanced by I/77 for 1 — 77 < t < 1. Taking if < 1 — rj in the bulk region, 
squaring the resulting differential operators, and adopting the Higgs profile given in (27), we 
thus need to solve 



d 2 


v 2 


dt 2 




d 2 


v 2 


dt 2 


2MI kV 2 



q q 






vS]y c 

q q 



q q 







yCfyS 

q q 



Al L (t,t';-p 2 )=0 



(C.6) 



where the dots denote subleading terms. The solutions to these equations involve hyperbolic 
trigonometric functions, whose arguments contain the matrices 



V2M } 



KK 



V2M } 



KK 



(C.7) 



and their hermitian conjugates. It is then not difficult to show that, in the limit 77 — > 0, the 
boundary conditions given in (57) still hold, provided we use X q as defined here instead of 
the original definition in (30), and Y q as shown in (72) instead of the original definition in 
(58). Solving the bulk equations of motion for the propagator functions with these boundary 
conditions, we recover our previous solutions with the substitutions just described. 



50 



References 

[1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 1 (2012) [arXiv: 1207.7214 [hep- 
ex]]. 

[2] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012) [arXiv:1207.7235 
[hep-ex]]. 

[3] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) [hep-ph/9905221]. 

[4] Y. Grossman and M. Neubert, Phys. Lett. B 474, 361 (2000) [hep-ph/9912408]. 

[5] T. Gherghetta and A. Pomarol, Nucl. Phys. B 586, 141 (2000) [hep-ph/0003129]. 

[6] S. J. Huber and Q. Shafi, Phys. Lett. B 498, 256 (2001) [hep-ph/0010195]. 

[7] K. Agashe, G. Perez and A. Soni, Phys. Rev. Lett. 93, 201804 (2004) [hep-ph/0406101]. 

[8] K. Agashe, G. Perez and A. Soni, Phys. Rev. D 71, 016002 (2005) [hep-ph/0408134]. 

[9] C. Csaki, A. Falkowski and A. Weiler, JHEP 0809, 008 (2008) [arXiv:0804.1954 [hep-ph]]. 

[10] S. Casagrande, F. Goertz, U. Haisch, M. Neubert and T. Pfoh, JHEP 0810, 094 (2008) 
[arXiv:0807.4937 [hep-ph]]. 

[11] M. Blanke, A. J. Buras, B. Duling, S. Gori and A. Weiler, JHEP 0903, 001 (2009) 
[arXiv:0809.1073 [hep-ph]]. 

[12] M. Blanke, A. J. Buras, B. Duling, K. Gemmler and S. Gori, JHEP 0903, 108 (2009) 
[arXiv:0812.3803 [hep-ph]]. 

[13] M. Bauer, S. Casagrande, U. Haisch and M. Neubert, JHEP 1009, 017 (2010) 
[arXiv:0912.1625 [hep-ph]]. 

[14] A. Djouadi and G. Moreau, Phys. Lett. B 660, 67 (2008) [arXiv:0707.3800 [hep-ph]]. 

[15] A. Falkowski, Phys. Rev. D 77, 055018 (2008) [arXiv:0711.0828 [hep-ph]]. 

[16] G. Cacciapaglia, A. Deandrea and J. Llodra-Perez, JHEP 0906, 054 (2009) 
[arXiv:0901.0927 [hep-ph]]. 

[17] G. Bhattacharyya and T. S. Ray, Phys. Lett. B 675, 222 (2009) [arXiv:0902.1893 [hep- 
ph]]. 

[18] C. Bouchart and G. Moreau, Phys. Rev. D 80, 095022 (2009) [arXiv:0909.4812 [hep-ph]]. 

[19] S. Casagrande, F. Goertz, U. Haisch, M. Neubert and T. Pfoh, JHEP 1009, 014 (2010) 
[arXiv: 1005.4315 [hep-ph]]. 

[20] A. Azatov, M. Toharia and L. Zhu, Phys. Rev. D 82, 056004 (2010) [arXiv:1006.5939 
[hep-ph]]. 

51 



[21] A. Azatov and J. Galloway, Phys. Rev. D 85, 055013 (2012) [arXiv:1110.5646 [hep-ph]]. 

[22] F. Goertz, U. Haisch and M. Neubert, Phys. Lett. B 713, 23 (2012) [arXiv:1112.5099 
[hep-ph]]. 

[23] M. Carena, S. Casagrande, F. Goertz, U. Haisch and M. Neubert, JHEP 1208, 156 (2012) 
[arXiv: 1204.0008 [hep-ph]]. 

[24] A. Azatov, M. Toharia and L. Zhu, Phys. Rev. D 80, 035016 (2009) [arXiv:0906.1990 
[hep-ph]]. 

[25] C. Delaunay, J. F. Kamenik, G. Perez and L. Randall, JHEP 1301, 027 (2013) 
[arXiv: 1207.0474 [hep-ph]]. 

[26] K. Agashe, A. Delgado, M. J. May and R. Sundrum, JHEP 0308, 050 (2003) [hep- 
ph/0308036]. 

[27] C. Csaki, C. Grojean, L. Pilo and J. Terning, Phys. Rev. Lett. 92, 101802 (2004) [hep- 
ph/0308038]. 

[28] K. Agashe, R. Contino, L. Da Rold and A. Pomarol, Phys. Lett. B 641, 62 (2006) [hep- 
ph/0605341]. 

[29] M. S. Carena, E. Ponton, J. Santiago and C. E. M. Wagner, Nucl. Phys. B 759, 202 

(2006) [hep-ph/0607106]. 

[30] G. Cacciapaglia, C. Csaki, G. Marandella and J. Terning, Phys. Rev. D 75, 015003 (2007) 
[hep-ph/0607146]. 

[31] R. Contino, L. Da Rold and A. Pomarol, Phys. Rev. D 75, 055014 (2007) [hep- 
ph/0612048]. 

[32] M. S. Carena, E. Ponton, J. Santiago and C. E. M. Wagner, Phys. Rev. D 76, 035006 

(2007) [hep-ph/0701055]. 

[33] R. Contino, Y. Nomura and A. Pomarol, Nucl. Phys. B 671, 148 (2003) [hep-ph/0306259]. 

[34] K. Agashe, R. Contino and A. Pomarol, Nucl. Phys. B 719, 165 (2005) [hep-ph/0412089]. 

[35] B. Batell, T. Gherghetta and D. Sword, Phys. Rev. D 78, 116011 (2008) [arXiv:0808.3977 
[hep-ph]]. 

[36] J. A. Cabrer, G. von Gersdorff and M. Quiros, New J. Phys. 12, 075012 (2010) 
[arXiv:0907.5361 [hep-ph]]. 

[37] J. A. Cabrer, G. von Gersdorff and M. Quiros, Phys. Lett. B 697, 208 (2011) 
[arXiv: 101 1.2205 [hep-ph]]. 

[38] P. R. Archer, JHEP 1209, 095 (2012) [arXiv: 1204.4730 [hep-ph]]. 



52 



[39] L. Randall and M. D. Schwartz, JHEP 0111, 003 (2001) [hep-th/0108114]. 

[40] A. Pomarol, Phys. Rev. Lett. 85, 4004 (2000) [hep-ph/0005293]. 

[41] K. -w. Choi, H. D. Kim and I. -W. Kim, JHEP 0211, 033 (2002) [hep-ph/0202257]. 

[42] W. D. Goldberger and I. Z. Rothstein, Phys. Rev. Lett. 89, 131601 (2002) [hep- 
th/0204160]. 

[43] K. Agashe, A. Delgado and R. Sundrum, Nucl. Phys. B 643, 172 (2002) [hep-ph/0206099]. 

[44] P. R. Archer, M. Carena, A. Carmona and M. Neubert, in preparation. 

[45] M. Puchwein and Z. Kunszt, Annals Phys. 311, 288 (2004) [hep-th/0309069]. 

[46] R. Contino and A. Pomarol, JHEP 0411, 058 (2004) [hep-th/0406257]. 

[47] M. S. Carena, A. Delgado, E. Ponton, T. M. P. Tait and C. E. M. Wagner, Phys. Rev. D 
71, 015010 (2005) [hep-ph/0410344]. 

[48] C. Csaki, Y. Grossman, P. Tanedo and Y. Tsai, Phys. Rev. D 83, 073002 (2011) 
[arXiv: 1004.2037 [hep-ph]]. 

[49] S. J. Huber, Nucl. Phys. B 666, 269 (2003) [hep-ph/0303183]. 

[50] F. Goertz, arXiv: 11 12.6387 [hep-ph]. 

[51] M. Beneke and M. Neubert, Nucl. Phys. B 651, 225 (2003) [hep-ph/02 10085]. 
[52] A. Djouadi, Phys. Rept. 459, 1 (2008) [hep-ph/0503173]. 

[53] M. E. Albrecht, M. Blanke, A. J. Buras, B. Duling and K. Gemmler, JHEP 0909, 064 
(2009) [arXiv:0903.2415 [hep-ph]]. 

[54] A. D. Medina, N. R. Shah and C. E. M. Wagner, Phys. Rev. D 76, 095010 (2007) 
[arXiv:0706.1281 [hep-ph]]. 

[55] Talks presented by F. Hubaut (ATLAS Collaboration) and M. G. Gomez-Ceballos (CMS 
Collaboration) at the Rencontres de Moriond, Electroweak Interactions and Unified 
Theories, La Thuile, Aosta Valley (Italy), 2-9 March 2013; slides for download at: 
https : //indico . in2p3. f r/conf erenceOtherViews .py?view=standard&conf Id=7411 

[56] M. E. Peskin and T. Takeuchi, Phys. Rev. D 46, 381 (1992). 

[57] M. Baak, M. Goebel, J. Haller, A. Hoecker, D. Kennedy, R. Kogler, K. Moenig and 
M. Schott et al, Eur. Phys. J. C 72, 2205 (2012) [arXiv:1209.2716 [hep-ph]]. 

[58] M. S. Carena, A. Delgado, E. Ponton, T. M. P. Tait and C. E. M. Wagner, Phys. Rev. D 
68, 035010 (2003) [hep-ph/0305188]. 



53 



[59] H. Davoudiasl, G. Perez and A. Soni, Phys. Lett. B 665, 67 (2008) [arXiv:0802.0203 
[hep-ph]]. 

[60] M. Bauer, R. Malm and M. Neubert, Phys. Rev. Lett. 108, 081603 (2012) 
[arXiv:1110.0471 [hep-ph]]. 

[61] M. Reece, arXiv: 1208. 1765 [hep-ph]. 

[62] M. Frank, N. Pourtolami and M. Toharia, arXiv:1301.7692 [hep-ph]. 

[63] G. Cacciapaglia, C. Csaki, G. Marandella and J. Terning, JHEP 0702, 036 (2007) [hep- 
ph/0611358]. 

[64] P. Breitenlohner and D. Z. Freedman, Annals Phys. 144, 249 (1982). 

[65] H. Davoudiasl, J. L. Hewett and T. G. Rizzo, Phys. Lett. B 473, 43 (2000) [hep- 
ph/9911262]. 



54 



